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Abstract 


MCMULLAN,  RICHARD  JEFFREY.  Numerical  Boundary  Conditions  Simulating 
the  Interaction  Between  Upstream  Disturbances  and  an  Axial  Compressor.  (Under 
the  direction  of  Dr.  D.  Scott  McRae.) 

New  small  disturbance  and  area  reduction  compressor  face  boundary  conditions 
that  model  the  unsteady  interactions  of  acoustic  disturbances  with  an  axial  compres¬ 
sor  are  presented.  The  new  small  disturbance  boundary  condition  is  formulated  to 
correct  the  deficiencies  associated  with  the  Paynter  small  disturbance  model.  The 
area  reduction  boundary  condition  provides  a  simple  approach  for  the  inlet  outflow 
boundary.  Both  of  these  boundary  conditions  are  implemented  in  one-dimensional 
and  axisymmetric  turbulent  flow  models  of  the  inlet/compressor  experiment  at  the 
University  of  Cincinnati.  Acoustic  reflections  from  the  compressor  face  boundary  con¬ 
ditions  are  compared  against  the  measured  experimental  reflection  characteristics  of 
the  axial  compressor.  The  performance  of  the  boundary  conditions  is  also  compared 
against  existing  boundary  conditions  such  as  the  Paynter  small  disturbance  boundary 
condition.  These  comparisons  show  that  the  new  small  disturbance  boundary  condi¬ 
tion  provides  the  best  accuracy  in  terms  of  the  prediction  of  the  reflected  disturbance 
from  the  interaction  of  an  acoustic  disturbance  with  a  compressor.  The  results  also 
show  that  the  area  reduction  boundary  condition  produces  acoustic  reflections  that 
agree  well  with  the  experimental  data. 
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1  Introduction 

The  success  of  commercial  supersonic  transport  programs  such  as  the  NASA  High 
Speed  Civil  Transport  (HSCT)  program  is  dependent  on  the  development  of  a  high- 
performance  propulsion  system.  For  flight  Mach  numbers  above  2.0,  a  mixed-compression 
inlet  is  needed  to  provide  high  total  pressure  recovery  and  low  flow  distortion  for  the 
propulsion  system.  A  mixed-compression  inlet  uses  a  system  of  oblique  shocks  up¬ 
stream  of  the  throat  and  a  normal  shock  just  aft  the  throat  to  decrease  the  inlet  flow 
Mach  number  from  supersonic  at  the  freestream  to  subsonic  at  the  face  of  the  ax¬ 
ial  compressor.  Supersonic  mixed-compression  inlets  are  susceptible  to  disturbances 
generated  in  the  atmosphere  (such  as  pressure  and  temperature  changes).  These 
disturbances  can  reduce  the  throat  and  normal  shock  Mach  numbers  and  cause  the 
shock  system  to  be  expelled  from  the  inlet  if  the  throat  Mach  number  falls  below  1.0 
or  the  normal  shock  Mach  number  is  less  than  the  throat  Mach  number.  This  unde¬ 
sirable  process  is  known  as  inlet  unstart.  The  penalties  of  inlet  unstart  are  significant 
reduction  in  total  pressure  recovery,  dramatic  increase  in  drag,  and  potential  loss  of 
aircraft  control.fi] 

In  the  past,  almost  all  of  the  performance  data  needed  to  design  the  inlet  was 
gathered  from  wind-tunnel  testing.  However,  increased  costs  of  wind-tunnel  testing 
have  led  to  the  use  of  computational  fluid  dynamic  (CFD)  techniques  for  the  design  of 
inlets.  In  order  to  properly  predict  inlet  unstart,  unsteady  CFD  codes  that  correctly 
model  the  propagation  of  disturbances  through  the  inlet  are  needed. 

Several  researchers  have  numerically  modeled  the  mixed-compression  inlet  and 
compressor  together  as  a  system  to  investigate  the  unsteady  flow  field  in  the  inlet. 

In  most  of  these  studies,  the  compressor  was  modeled  in  a  simple  one-dimensional 
sense. [2,  3,  4]  Suresh,  et.  al.  coupled  three-dimensional  inlet  and  turbomachinery 
CFD  codes  together  to  model  the  inlet/compressor  combination. [5]  However,  they 
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concluded  that  significant  computer  speed  up  would  be  required  to  make  the  cou¬ 
pling  of  inlet  and  compressor  together  a  viable  numerical  design  and  analysis  tool. 
Therefore,  supersonic  mixed-compression  inlet  flows  are  investigated  separately  from 
the  compressor  with  the  compressor  dynamics  represented  by  an  outflow  boundary 
condition. 

Atmospheric  disturbances  can  be  decomposed  into  acoustic  and  convective  dis¬ 
turbances  to  the  inlet  flow.  [6]  Acoustic  disturbances  are  pressure  disturbances  that 
propagate  downstream  through  the  inlet  at  the  flow  speed  plus  the  speed  of  sound. 
Convective  disturbances  are  thermodynamic  changes  that  propagate  downstream 
through  the  inlet  at  the  flow  speed.  Both  types  of  disturbances  eventually  inter¬ 
act  with  the  axial  compressor  and  reflect  back  upstream  as  acoustic  disturbances. 
Therefore,  a  viable  outflow  boundary  condition  for  an  unsteady  inlet  CFD  code  must 
properly  model  the  wave  reflection  from  the  compressor. 

Since  the  supersonic  flow  in  the  inlet  is  decelerated  to  subsonic  flow  at  the  com¬ 
pressor,  the  outlet  boundary  condition  at  the  compressor  face,  for  the  inlet  considered 
alone,  is  for  subsonic  outflow.  In  that  case,  one  physical  condition  must  be  specified 
at  the  boundary  and  the  remaining  conditions  are  calculated  using  the  information 
from  the  interior  of  the  inlet.  Traditionally,  subsonic  outflow  boundary  conditions 
have  been  evaluated  by  setting  a  selected  variable  to  be  constant  on  the  boundary 
and  solving  the  positive  eigenvalue  compatibility  equations  at  the  boundary.  This 
technique  is  mathematically  well  posed  but  does  not  model  the  presence  of  the  com¬ 
pressor. 

A  number  of  boundary  conditions  that  simulate  the  presence  of  the  compressor 
have  been  formulated  by  holding  a  variable  constant  at  the  boundary.  The  most 
commonly  selected  variables  have  been  constant  static  pressure[4,  7],  constant  axial 
velocity[4],  constant  axial  Mach  number[l,  8],  constant  volumetric  flow  rate[7],  and 
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constant  corrected  mass  flow  rate[l,  7].  All  of  these  boundary  conditions  are  easy  to 
implement  and  work  well  for  steady  inlet  flows. 

The  experimental  work  by  Freund  and  Sajben  was  conducted  in  a  facility  that 
mated  a  constant  area  annular  duct  with  a  multi-stage  axial  compressor  from  a  Gen¬ 
eral  Electric  helicopter  engine. [9, 10, 11, 12, 16, 17]  A  downstream  travelling  expansion 
pulse  was  created  from  the  collapse  of  a  flexible  bump.  The  reflection  from  the  axial 
compressor  was  an  expansion  pulse  travelling  upstream.  Four  pressure  transducers 
were  used  to  record  the  passage  of  the  incident  wave  from  the  bump  collapse  and 
the  reflected  wave  from  the  compressor.  The  experiment  provides  comparative  data 
for  testing  numerical  compressor  face  boundary  conditions.  In  Section  2,  a  detailed 
description  of  the  experiment  and  the  results  observed  will  be  presented. 

When  compared  with  the  experimental  data  collected  by  Freund  and  Sajben,  a 
constant  static  pressure  boundary  condition  was  shown  to  reflect  the  expansion  wave 
travelling  through  the  inlet  with  the  wrong  sign. [9, 10, 11, 12]  Although  constant  axial 
velocity,  constant  axial  Mach  number,  constant  volumetric  flow  rate,  and  constant 
corrected  mass  flow  rate  all  provide  the  proper  sign  of  the  reflection  wave,  the  compu¬ 
tational  results  showed  that  these  boundary  conditions  overpredict  the  magnitude  of 
the  acoustic  reflection  when  compared  with  the  experimental  data.  Therefore,  these 
simple  boundary  conditions  do  not  accurately  model  the  unsteady  behavior  at  the 
compressor. 

Mayer  and  Paynter  developed  a  boundary  condition  model  that  specified  com¬ 
pressor  face  corrected  mass  flow  rate  as  a  linear  function  of  the  stagnation  conditions.  [13] 
This  boundary  condition  was  found  to  provide  better  inlet  unstart  predictions  than 
setting  the  corrected  mass  flow  rate  as  a  constant.  In  addition,  Chung  and  Cole 
formulated  a  new  boundary  condition  that  specified  a  uniform  Mach  number  at  the 
compressor  face  and  allowed  the  static  pressure  to  vary  along  the  face. [14]  Numerical 
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calculations  performed  by  Slater  and  Paynter  showed  that  both  of  these  boundary 
conditions  grossly  overpredict  the  magnitude  of  the  reflected  acoustic  wave  when 
compared  to  University  of  Cincinnati  experiment. [15]  Therefore,  these  boundary  con¬ 
ditions  do  not  properly  model  the  unsteady  flow  reflection  at  the  compressor. 

Paynter,  et.  al.  have  developed  a  small  disturbance  model  for  the  interaction  of 
freestream  disturbances  with  an  axial  compressor. [18,  19,  20]  This  model  provides 
the  best  agreement  so  far  with  the  University  of  Cincinnati  experimental  data.  [15, 
20]  However,  this  model  does  not  take  into  account  the  turning  of  the  flow  by  a 
compressor’s  inlet  guide  vanes.  In  addition,  this  small  disturbance  boundary  condition 
assumes  that  no  work  is  being  done  by  the  compressor  at  steady  flow  conditions.  The 
Paynter  model  formulation  and  problems  will  also  be  presented  in  Section  2. 

The  current  work  uses  the  small  disturbance  ideas  of  Paynter,  et.  al.  to  develop 
a  new  small  disturbance  boundary  condition  that  accounts  for  the  presence  of  inlet 
guide  vanes  and  models  more  closely  the  compressor  dynamics.  Section  3  presents  the 
development  of  this  new  compressor  face  boundary  condition.  In  addition,  a  simple 
area  reduction  boundary  condition  was  developed  for  the  compressor  face  outflow 
boundary.  The  area  reduction  condition  compares  the  open  frontal  area  a  wave  sees 
as  it  enters  the  blade  passage  to  the  physical  cross-sectional  area  of  the  blade  passage. 
The  development  of  this  boundary  condition  is  also  presented  in  Section  3. 

Section  4  presents  the  governing  equations  for  both  the  one-dimensional  and  ax- 
isymmetric  two-dimensional  models  used  to  simulate  numerically  the  experimental 
work  performed  at  the  University  of  Cincinnati.  In  Section  5,  the  numerical  algo¬ 
rithms  used  to  solve  the  two  models  are  described  in  detail.  Section  6  compares  the 
performance  of  the  new  small  disturbance  boundary  condition  and  area  reduction 
boundary  condition  to  the  experimental  data  from  the  University  of  Cincinnati  and 
the  Paynter  small  disturbance  model  for  both  the  one-dimensional  and  axisymmetric 
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two-dimensional  models  of  the  experiment.  Sections  7  presents  the  conclusions  made 
about  the  use  of  these  new  boundary  conditions  in  unsteady  inlet  analysis. 
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2  Prior  Research 

2.1  Compressor  Face  Boundary  Condition  Experiment 

Freund  and  Sajben  at  the  University  of  Cincinnati  developed  an  experiment  to  study 
rapid  flow  transients  in  a  constant  area  annular  duct  mated  with  an  operating  axial 
compressor.  [9, 10, 11, 12, 16, 17]  The  experiment  used  a  collapsing  bump  to  create  one 
millisecond  duration  acoustic  expansion  pulses  travelling  both  upstream  and  down¬ 
stream  through  the  annular  duct.  The  amplitude  of  the  expansion  pulse  travelling 
downstream  toward  the  compressor  was  approximately  4%  of  the  mean  static  pres¬ 
sure.  The  incident  and  reflected  expansion  pulses  were  tracked  using  fast-response 
pressure  transducers.  From  the  pressure  transducer  data,  the  amplitude  and  dura¬ 
tion  of  the  reflected  pulse  were  determined.  This  experimental  research  provides  the 
CFD  community  with  a  database  to  validate  new  compressor  face  boundary  condition 
formulations.  However,  since  none  of  the  currently  available  inlet  outflow  boundary 
conditions  predict  correctly  the  data  collected  in  their  research,  Freund  and  Sajben 
concluded  that  new  outflow  boundary  conditions  are  needed.  [11,  12] 

2.1.1  Experimental  Apparatus  and  Test  Procedure 

The  experiment  at  the  University  of  Cincinnati  mated  a  constant  area  annular  duct  to 
an  axial  compressor.  Figure  1  shows  the  schematic  of  the  experiment  used  by  Freund 
and  Sajben  to  investigate  the  reflection  characteristics  of  the  axial  compressor  to  a 
downstream  travelling  acoustic  pulse. 

A  General  Electric  T58-3  helicopter  engine  with  the  power  turbine  removed  was 
used  for  this  experiment.  The  engine  was  modified  for  cold  operation  by  removing 
the  combustor  and  fuel  injection  systems  and  replacing  them  with  an  external  high 
pressure  air  supply  for  the  turbine.  The  maximum  attainable  speed  for  the  compressor 
in  cold  operation  was  approximately  20,000  RPM,  which  corresponds  to  an  axial 
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Modified  GET58  helicopter  engine 
Max.  speed  2  20,000  RPM 
Engine-face  Mach  number  a  0J20 


Spear  mechanism  g 
Vanes  m 
Structure  (she9)  ® 


Figure  1:  Inlet  Duct  and  Compressor  Rig  Schematic  (Dimensions  in  cm)  [11] 

compressor  face  Mach  number  of  approximately  0.2.  The  compressor  had  ten  stages 
and  a  pressure  ratio  up  to  2.7  for  the  laboratory  configuration. 

The  blade  configuration  data  for  the  compressor  is  considered  proprietary  in¬ 
formation  by  the  General  Electric  Company.  The  Department  of  Mechanical  and 
Aerospace  Engineering  at  North  Carolina  State  University  was  not  able  to  get  release 
of  this  information.  However,  the  following  descriptions  of  the  compressor  configura¬ 
tion  have  already  been  published. [11,  12,  15] 

The  compressor  had  ten  stages  and  a  variable  inlet  guide  vane  (VIGV)  row 
powered  by  a  linear  actuator.  The  VIGV  row  contained  38  uncambered  blades  whose 
stagger  angle  varied  from  31.7°  to  -5.3°.  The  first  stage  rotor  contained  30  blades 
with  a  solidity  and  stagger  angle  at  the  mean  radius  of  1.16  and  52°,  respectively.  In 
addition,  the  cross-sectional  area  was  reduced  across  the  VIGVs  by  87%. 

Figure  1  shows  the  1.80  m  long  constant  area  annular  duct  mated  with  the 
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axial  compressor.  The  duct’s  hub  and  case  diameters  were  0.137  m  and  0.258  m, 
respectively.  This  corresponded  to  an  annular  height  of  0.0605  m.  The  annular  duct 
had  a  flexible  collapsing  bump  whose  center  was  located  1.14  m  from  the  leading 
edge  of  the  VIGVs.  Figures  2  and  3  show  the  flexible  bump  deflated  and  inflated, 
respectively.  These  figures  also  show  that  the  transition  from  the  hub  to  the  bump 
for  both  cases  was  continuous.  The  bump  was  inflated  to  a  height  of  0.0111  m  at  its 
center  and  completely  deflated  in  0.0008  seconds  by  rupturing  a  diaphragm. 


Figure  2:  Deflated  Flexible  Bump  in  the  Constant  Area  Annular  Duct  [11] 


Figure  3:  Inflated  Flexible  Bump  in  the  Constant  Area  Annular  Duct  [11] 


Six  fast  response  pressure  transducers  were  used  to  monitor  the  travel  of  the 
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downstream  moving  incident  expansion  pulse  and  its  reflected  upstream  moving  ex¬ 
pansion  pulse.  Four  transducers  were  mounted  on  the  case  and  two  on  the  hub. 
Figure  4  shows  the  four  axial  locations  of  the  transducers. 


Figure  4:  Location  of  the  Pressure  Transducers  (Dimensions  in  m)  [11] 

The  following  test  procedure  was  employed  by  Freund  and  Sajben  for  the  exper¬ 
iments.  First,  the  flexible  bump  was  inflated  to  the  proper  height  using  helium  gas. 
Once  steady  state  conditions  were  achieved  in  the  inlet,  the  steady  state  flow  param¬ 
eters  of  ambient  pressure  and  temperature,  corrected  mass  flow  rate,  and  corrected 
compressor  speed  were  recorded.  Then,  the  diaphragm  was  burst  which  simulta¬ 
neously  caused  the  bump  to  deflate  and  the  pressure  transducers  to  acquire  data. 
Finally,  the  data  from  the  pressure  transducers  was  reduced  to  show  the  travel  of  the 
acoustic  waves  in  the  inlet. 

2.1.2  Experimental  Results 

Freund  and  Sajben  performed  many  experimental  tests  using  the  inlet  and  axial  com¬ 
pressor  set-up  at  low  speed  inlet  Mach  numbers  of  0.10  to  0.18.  The  data  collected 
from  these  runs  are  available  at  the  University  of  Cincinnati  Aerospace  Engineer¬ 
ing  anonymous  ftp  server  (ftp. ase.uc.edu). [11]  The  reduced  data  files  provided  the 
pressure  perturbation  histories  of  all  six  pressure  transducers  for  each  experimental 


test. 
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Figure  5  shows  the  plots  of  differential  pressures  at  each  of  the  four  transducer 
locations  on  the  annular  case  versus  time  for  an  inlet  Mach  number  Miniet  of  0.17. 
The  differential  pressures  5p  were  the  differences  in  pressure  from  the  steady  state 
pressure  recorded  at  each  transducer.  The  downstream  travelling  incident  expansion 
pulse  created  by  the  deflated  bump  is  labelled  IEP.  The  reflection  of  this  incident 
pulse  is  an  upstream  travelling  expansion  pulse  labelled  REP.  In  addition,  the  de¬ 
flated  bump  creates  an  upstream  travelling  expansion  pulse  which  is  reflected  off 
the  bellmouth  entrance  as  a  compression  pulse  labelled  RCP.  At  pressure  transducer 
4,  the  incident  and  reflected  expansion  pulses  overlap.  In  order  to  investigate  the 
shape  and  amplitude  of  the  reflected  expansion  pulse,  the  incident  pulse  measured  at 
transducer  1  was  shifted  in  time  and  subtracted  from  the  pressure  data  collected  at 
transducer  4. 

Figure  6  shows  an  expanded  view  of  the  reflected  expansion  pulse  in  Figure  5  at 
transducer  4.  The  maximum  amplitude  of  the  reflected  pulse  from  the  compressor 
is  approximately  36%  of  the  incident  pulse  amplitude.  From  plots  similar  to  this 
figure,  Freund  and  Sajben  determined  that  the  first  and  largest  negative  peak  in  the 
pulse’s  pressure  plot  represents  the  reflection  of  the  incident  pulse  from  the  first  stage 
rotor.  [11]  They  also  deduced  that  the  other  negative  peaks  were  reflections  from  the 
next  few  stages  in  the  compressor.  In  addition,  numerical  work  performed  by  Slater, 
Freund,  and  Sajben  showed  that  the  traditional  subsonic  outflow  boundary  conditions 
provided  reflection  pulses  that  differed  sharply  from  those  observed  from  experimental 
data,  for  example  the  reflected  pulse  given  in  Figure  6. [11,  17]  This  result  emphasizes 
the  need  for  the  development  of  new  compressor  face  boundary  conditions  that  agree 
more  closely  with  the  experimental  data. 

After  investigating  the  experimental  data,  Freund  and  Sajben  also  concluded  that 
the  axial  Mach  number  at  the  compressor  face  and  the  stagger  angle  of  the  VIGVs 


Pressure  Transducer  1 


Pressure  Transducer  2 
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Figure  6:  Expanded  View  of  Reflection  Pulse  at  Trnasducer  4  for  Miniet  -  0.17 

had  a  significant  effect  on  the  reflection  process  of  the  incident  expansion  pulse.  [11] 
The  amplitude  of  the  reflected  wave  increased  for  both  increases  in  the  Mach  number 
and  the  stagger  angle. 

2.2  Paynter  Small  Disturbance  Compressor  Face  Boundary 
Condition 

Paynter  et.  al.  have  developed  a  new  compressor  face  boundary  condition  based  on  a 
small  disturbance  model  of  the  unsteady  flow  around  a  compressor’s  first  stage.  [15, 
19,  20]  The  results  of  this  small  disturbance  boundary  condition  were  compared  to 
the  experimental  data  collected  at  the  University  of  Cincinnati.  Figure  7  shows  that 
this  compressor  face  boundary  condition  provides  significant  improvement  over  the 
existing  boundary  conditions  of  constant  pressure,  velocity,  and  Mach  number. 

2.2.1  Paynter  Small  Disturbance  Formulation 

Paynter  conducted  a  parametric  numerical  study  on  a  two-dimensional  cascade  to 
determine  the  amplitude  of  the  reflected  disturbance  as  a  function  of  the  blade  cam- 
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Figure  7:  Comparison  of  Numerical  and  Experimental  Results[20] 

ber,  blade  loading,  blade  solidity,  blade  stagger  angle,  incident  acoustic  disturbance 
strength,  and  upstream  axial  Mach  number. [18]  His  results  showed  that  only  the 
blade  stagger  angle  and  the  upstream  axial  Mach  number  had  a  strong  effect  on  the 
reflected  wave. 

Using  the  results  from  the  parametric  cascade  study,  Paynter  et.  al.  developed  a 
small  disturbance  boundary  condition  based  on  the  following  assumptions.  [20] 

1.  Blade  geometry  was  modeled  as  a  flat  plate  with  zero  thickness  and  a  solidity 
greater  than  1.0. 

2.  Flow  properties  prior  to  the  incident  acoustic  disturbance  were  known. 

3.  Prior  to  the  acoustic  disturbance,  the  flow  was  aligned  with  the  first  stage  rotor 
such  that  the  flow  properties  upstream  and  downstream  of  the  blade  passage 
were  the  same.  Figure  8  shows  a  schematic  of  the  flow  field  prior  to  the  inter¬ 
action  of  the  disturbance  with  the  rotor. 
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4.  Incident  acoustic  disturbance  was  a  step  change  in  static  pressure  and  travelled 
downstream  towards  the  blade  passage. 

5.  The  response  of  the  blade  passage  to  the  interaction  with  the  incident  acous¬ 
tic  disturbance  was  a  reflected  acoustic  disturbance  travelling  upstream  and 
a  transmitted  acoustic  disturbance  travelling  downstream.  Figure  9  shows  a 
schematic  of  the  flow  field  around  the  first  stage  rotor  after  its  interaction  with 
the  incident  disturbance. 

6.  Flow  at  the  exit  of  the  blade  passage  was  aligned  both  before  and  after  the 
transmitted  acoustic  disturbance. 


Figure  8:  Flow  Schematic  Prior  to  Acoustic  Disturbance [20] 


Figure  9:  Flow  Schematic  After  Interaction  of  Disturbance  with  Rotor[20] 

The  formulation  of  the  Paynter  small  disturbance  boundary  condition  used  the 
following  small  perturbation  relationships  derived  by  Shapiro.  [21]  For  right-running 
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(downstream  travelling)  acoustic  waves, 

Sp  6u  8u  .  . 

—  =7 — =7  M —  (1) 

p  a  u 

For  left-running  (upstream  travelling)  acoustic  waves, 

tip  &U  /r>\ 

—  =  -7 —  =  -7  M —  (2) 

pa  u 

In  addition,  Paynter  et.  al.  used  the  following  notation  system  for  relating  the 
flow  properties  in  regions  1  —  4  of  Figures  8  and  9. 

P2=P\JrP,2  u2  =  u1+u'2 

P3=Pi+  Pz  u3  =  Ui  +  u'z  (3) 

Pa  =  Pi  +  Pa  ua  =  U]_  +  u\ 


Region  1  was  designated  as  the  reference  state,  and  the  prime  quantities  represented 
small  perturbations  from  the  reference  state. 

Using  Equations  (2)  and  (3),  the  flow  properties  across  the  reflected  acoustic 
disturbance  were  given  by 


Pz 


Uo 


—  +  7-^xi  • 

Pi  Ul 


2— 

Pi 


(4) 


Assuming  no  time  rate  of  change  of  mass  in  a  control  volume  fixed  to  the  blade 
passage  given  in  Figure  9,  the  continuity  equation  for  the  control  volume  was  given 
by 

?3+7!4_^_7<=0  (5) 

Pi  Ul  Pi  u  1 

Using  this  control  volume  and  assuming  no  time  rate  of  change  of  energy  within 
the  control  volume,  the  energy  equation  was  given  by 


^3+^M2  =Q 

Pi  XlU\  Pi  7cos2rwi 


(6) 
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Finally,  the  flow  properties  across  the  transmitted  acoustic  disturbance  were 


given  by 

PA-^<=  0  (7) 

Pi  COS  r  Ml 

Equations  (4)  through  (7)  formed  a  system  of  linear  equations  which  were  solved 
to  determine  a  formulation  for  the  reflection  response  coefficient.  This  response  co¬ 
efficient  /?  was  defined  to  be  the  ratio  of  reflected  acoustic  disturbance  strength  to 


incident  acoustic  disturbance  strength. 

^=^P2=tan2 


6) 


1  +  Mx  i 
1  —  Mx  i 


.  (8) 

P2-P1 

Paynter  et.  al.  used  the  linearized  Euler  equations  in  characteristic  form  to  apply 
the  small  disturbance  boundary  condition  into  a  CFD  code.  [20]  Using  region  1  in 


Figures  8  and  9  as  the  reference  state,  the  linearized  Euler  equations  in  characteristic 
form  were  given  by 
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After  investigating  the  characteristic  behavior  of  the  boundary  where  the  incident 


acoustic  disturbance  interacts  with  the  compressor,  Paynter  et.  al.  concluded  that 
r  +  (3q  was  a  Reimann  invariant  on  the  boundary.  [20]  Then,  the  small  disturbance 


compressor  face  boundary  equation  for  pressure  was  given  by 

^  =  -(«,  + o,)  (1  +  P)%  (10) 

where  /3  was  calculated  using  Equation  (8).  In  addition,  the  density  and  velocity  at 
the  boundary  were  calculated  using  the  positive  eigenvalue  characteristic  equations 
in  Equation  (9)  and  the  new  pressure  calculated  from  Equation  (10). 
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2.2.2  Problems  with  Paynter  Small  Disturbance  Formulation 

Paynter  et.  al.  applied  their  small  disturbance  boundary  condition  across  the  blade 
passage  of  a  compressor’s  first  stage  rotor.  The  only  two  parameters  that  affected 
their  reflection  response  coefficient  were  the  axial  Mach  number  before  the  incident 
disturbance  and  the  stagger  angle  of  the  rotor  blades.  They  assume  that  the  flow 
entering  the  rotor  was  purely  axial.  This  assumption  neglected  the  flow  turning 
provided  by  inlet  guide  vanes. 

In  the  Paynter  model,  the  rotor  blades  were  modeled  as  flat  plates  so  that  the 
initial  steady  conditions  on  both  sides  of  the  blade  passage  were  equal.  This  model 
of  the  rotor  did  zero  work  on  the  fluid  at  steady  state.  In  comparison,  rotor  blades 
in  a  real  compressor  have  camber,  and  the  first  stage  rotor  does  work  on  the  fluid.  In 
addition,  the  flat  plate  model  did  negative  work  on  the  fluid  (i.e.  acted  as  a  turbine) 
when  the  incident  acoustic  disturbance  was  a  compression  wave.  Therefore,  the  flat 
blades  inaccurately  model  actual  rotor  blades.  This  small  disturbance  model  needs 
to  be  corrected  to  consider  inlet  guide  vanes  and  real  blade  shapes. 

2.2.3  Modification  to  Paynter  Small  Disturbance  Formulation 

Sajben  developed  a  small  disturbance  compressor  face  boundary  condition  that  was 
very  similar  to  the  Paynter  model. [22]  The  only  exception  was  that  Sajben  assumed 
the  transmitted  disturbance  did  not  travel  along  the  direction  of  the  blade  passage. 
In  his  model,  the  transmitted  wave  propagated  along  the  axial  direction,  and  an 
additional  vorticity  wave  was  created  off  the  trailing  edge  of  the  blades  to  align  the  flow 
behind  the  transmitted  disturbance  with  the  blade  passage.  According  to  Sajben,  the 
change  in  flow  angle  into  the  blade  passage  created  lift  on  the  blades  which  produced 
a  starting  vortex  from  each  blade.  Using  this  new  assumption,  Sajben’s  formulation 
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for  the  reflection  response  coefficient  was  given  by 

8= _ M«'tanT _  (11) 

P  2  (1  +  Mxl)  +  Mxl  tan2  T  V  ’ 

Numerical  calculations  performed  by  Slater  and  Paynter  showed  that  the  Paynter 
small  disturbance  boundary  condition  predicted  more  closely  the  magnitude  of  the 
reflected  acoustic  disturbance  than  the  Sajben  model  when  compared  to  University 
of  Cincinnati  experimental  data.  [15]  The  results  in  Section  6  will  also  show  that  the 
Paynter  model  performs  better  than  the  Sajben  model.  In  addition,  the  Sajben  small 
disturbance  boundary  condition  modification  did  not  improve  on  any  of  the  problems 
that  occurred  with  the  Paynter  model. 
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3  New  Compressor  Face  Boundary  Conditions 

Two  new  compressor  face  boundary  conditions  were  developed  to  model  the  axial 
compressor  reflection  process.  These  boundary  conditions  must  accurately  calculate 
the  amplitude  of  the  reflected  disturbance  and  be  simple  enough  to  implement  into  ex¬ 
isting  CFD  codes.  In  this  section,  both  the  new  small  disturbance  and  area  reduction 
boundary  conditions  will  be  presented. 

3.1  New  Small  Disturbance  Boundary  Condition 

A  new  small  disturbance  compressor  face  boundary  condition  that  models  the  inter¬ 
action  of  an  acoustic  disturbance  and  an  axial  compressor  was  developed  to  correct 
the  problems  in  the  Paynter  small  disturbance  formulation.  This  model  accounts  for 
the  turning  of  the  flow  through  the  inlet  guide  vanes.  In  addition,  the  model  uses 
cambered  blade  geometry  and  the  velocity  of  the  compressor  blades  in  the  calcula¬ 
tion  of  the  reflection  response  coefficient  /?  from  an  axial  compressor.  Therefore,  this 
formulation  includes  more  of  an  axial  compressor’s  flow  dynamics. 

3.1.1  Small  Disturbance  Formulation 

Most  of  the  same  assumptions  employed  by  Paynter  et.  al.  were  also  used  in  this  new 
boundary  condition  development.  However,  changes  were  made  to  model  the  effects 
of  the  inlet  guide  vanes  and  the  cambered  blade  designs.  This  small  disturbance 
boundary  condition  is  formulated  based  on  the  following  assumptions: 

1.  The  flow  properties  prior  to  the  incident  acoustic  disturbance  interacting  with 
the  axial  compressor  are  known. 

2.  Prior  to  the  acoustic  disturbance,  the  flow  is  aligned  with  the  blade  passages  of 
both  the  inlet  guide  vanes  and  the  first  stage  rotor. 
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3.  The  flow  is  isentropic  through  the  inlet  guide  vanes. 

4.  The  incident  acoustic  disturbance  is  a  step  change  in  static  pressure  and  trav¬ 
elled  downstream  towards  the  compressor. 

5.  The  response  of  each  blade  passage  to  the  interaction  with  an  acoustic  distur¬ 
bance  is  a  reflected  acoustic  disturbance  travelling  upstream  and  a  transmitted 
acoustic  disturbance  travelling  downstream. 

6.  The  flow  at  the  exit  of  each  blade  passage  is  aligned  with  the  blade  angle  both 
before  and  after  its  interaction  with  the  acoustic  disturbance. 


This  boundary  condition  models  an  acoustic  disturbance’s  travel  through  the 
inlet  guide  vanes  and  first  stage  rotor  of  an  axial  compressor.  First,  the  acoustic 
disturbance  propagates  downstream  through  the  inlet  at  the  flow  speed  plus  the  speed 
of  sound  and  interacts  with  the  inlet  guide  vanes  of  the  axial  compressor.  Figure  10 
shows  schematically  the  flow  field  around  the  inlet  guide  vanes  prior  to  the  arrival  of 
the  incident  acoustic  disturbance.  The  flow  Mach  number  M\  behind  the  inlet  guide 
vanes  is  calculated  from  the  following  equation  for  constant  mass  flow  rate  using  a 
Newton  iteration. 


S\Mi  cos  o.i 


S0Mo 


1  + 


n  ^+1 
2(7-l) 


1  +  1— -M? 


7+1  . 
2(7-1) 


(12) 


where  S  is  the  cross-sectional  area  and  o.\  is  the  blade  passage  exit  angle.  Using 
the  assumption  of  isentropic  flow  across  the  inlet  guide  vanes,  the  thermodynamic 
variables  at  the  exit  of  the  guide  vanes  are  computed.  In  addition,  the  axial  Mach 
number  Mx i,  axial  velocity  iq,  and  circumferential  velocity  w\  are  calculated  at  the 
exit. 


The  result  of  the  interaction  of  the  acoustic  disturbance  with  the  inlet  guide  vanes 
is  a  reflected  wave  travelling  upstream  and  a  transmitted  wave  travelling  downstream 
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Acoustic 

Disturbance 


Figure  10:  Flow  Schematic  of  Inlet  Guide  Vanes  Prior  to  Incident  Acoustic  Distur¬ 
bance 

aligned  with  the  exit  angle  of  the  blade  passage.  Figure  11  depicts  the  flow  around 
the  inlet  guide  vanes  after  its  interaction  with  the  incident  acoustic  wave. 


Reflected 

Disturbance 


Transmitted 


Figure  11:  Flow  Schematic  of  Inlet  Guide  Vanes  After  Interaction  With  Incident 
Acoustic  Disturbance 


The  following  small  disturbance  notation  system  is  used  to  relate  the  flow  prop¬ 


erties  in  regions  0,  1,  2,  6,  and  7  of  Figures  10  and  11. 


Pe  -  Po  +  p'6  ue  =  u0  +  u'6  p6  =  p0+  Pq 

P7  —  P0  +  P7  u7  —  u0  +  u7  P7  =  Po  +  P7 

P2=Pl+  P2  U2  =  Ui+  u'2  p2=  pi  +  P2 


(13) 


Regions  0  and  1  are  designated  as  the  reference  states  in  front  of  and  behind  the  inlet 
guide  vanes,  respectively. 

Using  Equations  (1)  and  (13),  the  flow  properties  across  the  incident  acoustic 
disturbance  are  given  by 

*  =  (14) 

Po  Uo 

Using  Equations  (2)  and  (13),  the  flow  properties  across  the  reflected  acoustic 
disturbance  are  given  by 


a  -  a  -  -imm 


« 6 

\u0  u0 


(15) 


Po  Po 

This  equation  is  simplified  into  the  following  form  by  substituting  Equation  (14)  into 
it  and  rearranging. 

nrJ  nr! 

(16) 


+  jMxo—  =  2^6 


Po  u0  po 

Assuming  no  time  rate  of  change  of  mass  in  a  control  volume  attached  to  the 
inlet  guide  vanes,  the  continuity  equation  has  the  following  forms. 


P7U7S0  —  p2u2S\ 

P7U0S0  +  PqUjSq  =  p'2u\ S\  +  piu'2S\ 


(17) 


Using  the  isentropic  relation 


dp  1 8p 
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Equation  (17)  is  reduced  to  the  following  form. 


*4 +7?4  =  6.  +  7!f2 

Po  uo  Pi  ui 


(19) 


Assuming  no  time  rate  of  change  of  energy  for  this  control  volume,  the  conser¬ 
vation  of  energy  equation  is  given  by 


7  Pi 
7“  lp7 


h 


tl 


+ 


ht2 


7  P2 
7-  lp2 


+  -  («2  +  TO,) 


(20) 


where  ht  is  the  total  enthalpy.  Using  Equation  (18)  and  the  velocity  component 
relation  tu2  =  u2  tane^,  Equation  (20)  is  reduced  to  the  following  form. 


Pi  P2  ,  U2Ul 

- h  UjUp  — - 1 - 7. - 

PO  Pi  COS2  Q!i 


(21) 


Using  Equations  (1)  and  (13),  the  flow  properties  across  the  transmitted  acoustic 
disturbance  are  given  by 

nr!  r\i  1\/f  - 

(22) 


P2  _  lMxl  U2 
Pi 


COS  OKI  Wl 

Equations  (16),  (19),  (21),  and  (22)  form  a  system  of  linear  equations  for  four 
unknowns  p'7,  p'2,  u'r,  and  u'2.  This  linear  system  is  solved  to  determine  ratio  of  the 
transmitted  acoustic  disturbance  strength  to  the  known  incident  acoustic  disturbance 
strength. 

P2 

2a0coso!i(uo  +  «o) 


Pi  _ _ 

P1  ~ 

—  (ui  +  ai  cos  aU  (  eq  +  a0 
Po  1 


(;) 


COS  Ot\ 


(23) 


In  addition,  the  ratio  of  reflected  acoustic  disturbance  strength  to  the  known  incident 
disturbance  strength  is  given  by 


P?  _  Pe 

Po _ po_  _  (1  +  Mx0)(uiai  -  upOp  coso;!) 

Pe  (1  -  Mx0 ) (uiai  +  u0a0  cos  £*i) 
Po 


(24) 


24 


Comparing  with  the  data  from  the  University  of  Cincinnati  Mach  0.17  run  presented 
in  Figure  5,  the  reflected  disturbance  from  the  interaction  of  the  incident  disturbance 
with  the  inlet  guide  vanes  was  found  to  be  very  weak.  The  ratio  given  in  Equation 
(24)  was  calculated  to  be  approximately  0.02.  This  negligible  amount  of  reflection 
off  the  inlet  guide  vanes  agrees  with  the  University  of  Cincinnati  conclusion  that  the 
main  reflection  occurs  at  the  first  stage  rotor. 

After  transmission  through  the  inlet  guide  vanes,  the  acoustic  disturbance  travels 
downstream  towards  the  first  stage  rotor.  Figure  12  depicts  the  flow  around  the 
first  stage  rotor  prior  to  the  arrival  of  the  inlet  guide  vanes’  transmitted  acoustic 
disturbance.  The  relative  velocity  vector  measured  with  respect  to  the  moving  rotor 
vt  is  given  by 

\t  =  \?-l}  (25) 

where  is  the  absolute  velocity  measured  with  respect  to  the  fixed  reference  frame 
of  the  compressor  rig  and  ij  is  rotor’s  blade  velocity.  The  axial  component  of  the 
velocity  remains  the  same  in  both  reference  frames.  The  relative  flow  properties 
upstream  of  the  rotor  in  region  1  are  calculated  using  the  known  rotor  blade  speed  U 
and  the  previously  computed  absolute  flow  conditions. 

The  relative  Mach  number  M4r  downstream  of  the  first  stage  rotor  is  calculated 
using  the  constant  mass  flow  rate  equation  for  the  rotor  blade  passage. 


SiPnrMlr  cos  /?i 
1  - 

1  + 


1  r 


2(7-1) 


S'4pi4rM4r  cos  /?4 

!  _  -y+l 

7—1  „  1 2#Ay 

1  +  UAK 


(26) 


where  ptr  is  the  relative  total  pressure,  fix  is  blade  inlet  angle,  and  /?4  is  the  blade  exit 
angle.  The  relative  total  pressure  ratio  across  the  rotor  passage  is  given  by 

Pt4r 


Ptlr 


=  i-0, 


cr 


l Ml 

1  + 


1  r 


7-1 


(27) 
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Figure  12:  Flow  Schematic  of  First  Stage  Rotor  Prior  to  Arrival  of  Inlet  Guide  Vanes’ 
Transmitted  Acoustic  Disturbance 

where  (f)cr  is  the  rotor  loss  coefficient  with  a  typical  range  in  value  from  0.05  to  0.12. 
This  loss  coefficient  accounts  for  the  non-isentropic  flow  through  the  rotor  passage. 
Using  the  relative  total  pressure,  relative  Mach  number,  rotor  blade  speed,  and  exit 
blade  angle,  all  of  the  remaining  flow  variables  in  region  4  are  computed. 

The  result  of  the  interaction  of  the  transmitted  acoustic  disturbance  with  the 
first  stage  rotor  is  a  reflected  wave  travelling  upstream  aligned  with  the  exit  angle  of 
the  inlet  guide  vanes  and  a  transmitted  wave  travelling  downstream  aligned  with  the 
exit  angle  of  the  rotor.  Figure  13  depicts  the  flow  around  the  first  stage  rotor  after 
its  interaction  with  the  transmitted  acoustic  wave. 

The  following  small  disturbance  notation  system  is  used  to  relate  the  flow  prop¬ 
erties  in  regions  1,  2,  3,  4,  and  5  of  Figures  12  and  13. 

P2  =  Pi  +  P'2  U2  =  Ul+  u'2  p2  =  pi  +  p'2 

P3=Pl+  Pz  «3  =  Ui  +  u'z  p3=  pl+  p'3  (28) 

P5  =  Pa  +  p'5  u5  =  u4  +  u'5  p5  =  pA  +  p's 
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Figure  13:  Flow  Schematic  of  First  Stage  Rotor  After  Interaction  With  Inlet  Guide 
Vanes’  Transmitted  Acoustic  Disturbance 


Regions  1  and  4  are  designated  as  the  reference  states  in  front  of  and  behind  the  first 
stage  rotor,  respectively. 

Using  Equations  (2)  and  (28),  the  flow  properties  across  the  reflected  acoustic 
disturbance  from  the  rotor  are  given  by 


Ps 


P2 


lMxi 


/X 

\Ui 


«2 

Ml 


(29) 


Pl  Pi  COSCKi  \U  1 

This  equation  is  simplified  into  the  following  form  by  substituting  Equation  (22)  into 
it  and  rearranging. 

nrJ  m! 

(30) 


?3  ,  7A4l  U'3  _  2P2 

Pl  COS  Oil  U\  Pi 


Assuming  no  time  rate  of  change  of  mass,  the  continuity  equation  for  a  control 
volume  fixed  to  the  moving  first  stage  rotor  has  the  following  form. 


P3U1S1  +  piii  3S 1  —  p'3u/[Si  +  P4M5S4 


(31) 
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Using  Equation  (18)  on  either  side  of  the  rotor  passage,  Equation  (31)  is  reduced  to 
the  following  form. 


Ps 


Ur  po  Uo 

—  +  7—  =  —  +  7— 
P4  u4  pi  Ui 


(32) 


Assuming  no  time  rate  of  change  of  energy,  the  conservation  of  energy  equation 
for  this  control  volume  fixed  to  the  moving  rotor  is  given  by 


k 


thr 


=  h 


t3r 


_J_P 5 
7-lps 


+  o  (u5  +  wjr) 


7  P3 
7  -  lp3 


+  -(ul  +  wl) 


(33) 


where  w^r  =  u5  tan  j54  and  w3r  =  U  —  u3  tan  a4.  Using  Equation  (18)  on  both  sides  of 
the  rotor,  Equation  (33)  is  reduced  to  the  following  form. 


pk+<^  =  pA  +  j^_Uu, 

P4  COS2  P4  pi  COS2  ai 


(34) 


Using  Equations  (1)  and  (28),  the  flow  properties  across  the  transmitted  acoustic 
disturbance  from  the  rotor  are  given  by 


A  =  (35) 

p4  COS  P4  «4 

Equations  (30),  (32),  (34),  and  (35)  form  a  system  of  linear  equations  for  four 
new  unknowns  p'5,  p3,  u'5,  and  u'3.  This  linear  system  is  solved  to  determine  the 
ratio  of  the  rotor’s  reflected  acoustic  disturbance  strength  to  the  inlet  guide  vanes’ 
transmitted  acoustic  disturbance  strength. 


ti_P2 

Pi _ lh 

Pi 

Pi 


(cos  ai  +  Mxi)(u4a4  cos  a4  —  u\ai  cos  /34)  +  Uui  sin  a i  cos  ai  cos  /54 
(cos  ai  —  Mx\){u4a4  cos  +  «iOi  cos  /34)  +  Uu4  sin  a4  cos  a4  cos  p4 


(36) 


Equation  (36)  can  be  reduced  to  Equation  (8)  for  Paynter’s  reflection  response 
coefficient  if  the  turning  of  the  flow  by  the  inlet  guide  vanes  is  neglected  (0:4  =  0)  and 
the  blades  are  assumed  to  be  flat  plates  (uiaj  =  u4a4). 
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After  the  reflection  from  the  first  stage  rotor,  this  reflected  acoustic  disturbance 
travels  back  upstream  towards  the  inlet  guide  vanes.  Figure  14  shows  schematically 
the  flow  field  around  the  inlet  guide  vanes  prior  to  the  arrival  of  the  rotor’s  reflected 
acoustic  disturbance.  The  result  of  the  interaction  of  the  acoustic  disturbance  with 
the  inlet  guide  vanes  is  a  transmitted  wave  travelling  upstream  through  the  inlet 
and  a  reflected  wave  travelling  downstream  aligned  with  the  exit  angle  of  the  blade 
passage.  Figure  15  depicts  the  flow  around  the  inlet  guide  vanes  after  its  interaction 
with  the  rotor’s  reflected  acoustic  wave. 


Figure  14:  Flow  Schematic  of  Inlet  Guide  Vanes  Prior  to  Arrival  of  Rotor’s  Reflected 
Acoustic  Disturbance 

The  following  small  disturbance  notation  system  is  used  to  relate  the  flow  prop¬ 
erties  in  regions  2,  3,  7,  8,  and  9  of  Figures  14  and  15. 

P2  =  Pi  +  Pa  u2  =  tti  +  u'2  p2  =  pi  +  p'2 

PS=Pl  +  p's  U3  =  Ml  +  u'3  p3  =  px  +  p'z 

P9=Pl+  p'g  Ug  =  Ui+  u'9  p9=  pi+  p9  (37) 

Pi  =  Po  +  P?  u7  =  u0  +  u'7  p7  =  p0  +  p'? 

P8  =  PO  +  Ps  U&  =  U0+  u'8  p8  =  P0  +  P8 
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Figure  15:  Flow  Schematic  of  Inlet  Guide  Vanes  After  Interaction  With  Rotor’s 
Reflected  Acoustic  Disturbance 


Regions  0  and  1  are  still  designated  as  the  reference  states  in  front  of  and  behind  the 
inlet  guide  vanes,  respectively. 

Using  Equations  (1)  and  (37),  the  flow  properties  across  the  reflected  acoustic 
disturbance  are  given  by 


Pv  _  Ps  _  jMxi 

Pi  pi  cosax  \ui  Ui) 


(38) 


This  equation  is  simplified  into  the  following  form  by  substituting  Equation  (30)  into 
it  and  rearranging. 


Pg  TMsI  Mg 
Pi  COS  £*1  U\ 


=  2 


Ps 

Pi 


P_2 

Pi, 


(39) 


Using  the  assumption  of  no  time  rate  of  change  of  mass,  the  continuity  equation 
for  a  control  volume  fixed  to  the  inlet  guide  vanes  has  the  following  form. 


Psu  0*So  +  P0uS^0  —  P9U  1&1  +  Plu9^l  (40) 

Using  Equation  (18)  across  the  inlet  guide  vanes,  Equation  (40)  is  reduced  to  the 
following  form. 

p'o  u'o  p'a  u'a 
— +  7—  =  — +  7— 

Po  Uo  Pi  Ui 


(41) 
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Using  the  assumption  of  no  time  rate  of  change  of  energy,  the  conservation  of 
energy  equation  for  this  control  volume  fixed  to  the  blade  passage  is  given  by 


hts  —  htg 
7  P8  +  1  2  _  7  P9  1  (  2  ,  2) 

-lp«  +  2“8  "  7-lP9  2^9  +  ^ 


(42) 


7-ips  'i  ~  7  -  1  p9 

Using  Equation  (18)  and  the  velocity  component  relation  w$  =  u9  tan  cvi,  Equation 
(42)  is  reduced  to  the  following  form. 


Fs,  ,  _ P9  .  «9«1 

- b  U8Uo  — - 1 - 


(43) 

Po  ~  pi  COS*  ai 
Using  Equations  (2)  and  (37),  the  flow  properties  across  the  transmitted  acoustic 
disturbance  are  given  by 


*  -  a  = 

Po  Po 


_  «7 
\U0  U0 


(44) 


This  equation  is  simplified  into  the  following  form  by  substituting  Equation  (16)  into 
it  and  rearranging. 

i  =  (45) 

Po  U0  Po 

Equations  (39),  (41),  (43),  and  (44)  form  a  system  of  linear  equations  for  four  new 
unknowns  p'8,  p'(J,  u's,  and  u'9.  This  linear  system  is  solved  to  determine  an  equation 
for  the  strength  of  the  transmitted  acoustic  disturbance  travelling  upstream  through 
the  inlet. 


PO  _  [2(CQSQ!  +  Mxi)(u  1Q1  -  MpCOSCXi)]  Pe 

p0  [(1  -  Mx o)(cosai  +  Mxl)(u0a0  cos  ai  +  Uidi)]p0 
[2 a\Mxo  cos  ai  (cos  ai  -  M^) ] 

[(1  -  Mx 0)(cosa;i  +  Mxi)(u0a0  cosai  +  Uiai)] 


P0_P2_\ 
Pi  Pi) 


(46) 


Equations  (23)  and  (36)  are  used  to  simplify  the  previous  equation  into  an  equation 
for  the  compressor  response  coefficient  j3.  The  response  coefficient  for  this  new  small 
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disturbance  boundary  condition  is  given  by 


Ps-Pe 


^Miai  —  itoao  cosai 


+ 


1  -  Mx0 
1  +  Mx  o 


( 


Uidi  +  «oao  cos  Oi\ 

((cos ai  —  Mxl)  (4uju0«oOi cos2 c*i) 

(cosai  +  Mx i)2  (u\Cl\  +  u0°o  cos  ai)2 , 

(cos  cci  +  Mxi)(u4<i4  cos  —  U\0,\  cos  /?4)  +  Uu\  sin ot\  cos  a\  cos  /3 4 


(47) 


(cos  —  MIi)(lt404  COS  OL\  +  «iOi  COS  ^4)  +  UU\  sin  OL\  COS  Oil  COS  /?4 ; 

This  new  formulation  for  the  reflection  response  coefficient  takes  in  account  the  turn¬ 
ing  of  the  flow  through  the  inlet  guide  vanes  and  the  camber  and  speed  of  the  blades. 


3.1.2  Model  Implementation 

The  new  small  disturbance  compressor  face  boundary  condition  was  implemented 
into  a  CFD  code  similar  to  the  process  described  by  Paynter  et.  al.[20]  First,  the  new 
reflection  response  coefficient  f3  is  calculated  using  Equation  (47).  This  new  response 
coefficient  only  depends  on  the  flow  conditions  and  geometry  of  the  compressor  prior 
to  the  arrival  of  the  incident  acoustic  disturbance.  The  previous  time  level  conditions 
at  the  compressor  face  boundary  (region  0)  and  the  geometry  for  the  inlet  guide  vanes 
and  first  stage  rotor  varies  in  the  radial  direction.  Therefore,  the  needed  reference 
conditions  in  regions  1  and  4  of  Figures  10  and  12  and  [3  are  calculated  at  each 
computational  cell  along  the  outflow  boundary.  For  a  one-dimensional  model  of  the 
inlet,  mean  values  of  the  compressor’s  geometry  are  used  to  calculate  the  response 
coefficient. 

After  calculating  the  new  reflection  response  coefficient,  the  updated  pressure  at 
the  boundary  is  calculated  from  Equation  (10).  Then,  the  remaining  flow  variables 
at  the  compressor  face  boundary  are  updated  using  the  one  dimensional  positive 
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eigenvalue  compatibility  equations  given  by 


dp  1  dp  f  dp  1  dp\ 

dt  a2  dt  +  U  \  dx  a2  dx ) 

du  l  dp  .  .  ( du  1  dp\ 

dt  +  pa  dt  U  °  pa  dx  J 


=  0 
=  0 


and  the  equation  for  the  advection  of  the  radial  velocity  v  given  by 


dv  dv  n 

m  +“&=0 


3.2  Area  Reduction  Boundary  Condition 


(48) 


(49) 


For  this  simple  compressor  face  boundary  condition,  the  first  stage  rotor  blades  are 
modeled  as  flat  plates.  The  flow  in  front  of  the  downstream  travelling  acoustic  dis¬ 
turbance  is  assumed  to  be  aligned  with  the  rotor  blade  passage.  Figure  16  shows  the 
frontal  area  open  to  the  flow  if  the  flow  is  aligned  with  the  stagger  angle  T  of  the 
rotor  blades.  The  frontal  area  the  aligned  flow  saw  as  it  entered  the  blade  passage  is 
given  by 

ni  =  scosT  (50) 


Figure  16:  Frontal  Area  Open  to  Flow  Aligned  With  Blade  Passage 
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Behind  the  acoustic  disturbance,  the  flow  entering  the  blade  passage  is  not 
aligned  with  the  rotor  blades.  Therefore,  the  open  frontal  area  seen  by  the  flow 
is  reduced.  Figure  17  shows  the  reduced  open  frontal  area  for  a  flow  that  is  not 
aligned  with  the  blade  passage.  The  shaded  area  represents  the  frontal  area  that  the 
entering  flow  sees  as  closed.  The  frontal  area  that  is  closed  to  the  non-aligned  flow  is 
given  by 

n2  =  ctan  |/32  —  T|  (51) 

where  /?2  is  inlet  flow  angle. 


Figure  17:  Frontal  Area  Open  to  Flow  Not  Aligned  With  Blade  Passage 


The  area  reduction  boundary  condition  compares  the  open  frontal  area  the  flow 
after  the  disturbance  sees  to  that  seen  by  the  flow  in  front  of  the  acoustic  disturbance. 
Using  Equations  (50)  and  (51),  the  relative  amount  of  area  reduction  AR  is  given  by 


AR 


ni  -  n2 
ni 

^  ctan  |/32  -  T 
s  cosT 


(52) 


In  order  to  implement  this  area  reduction  boundary  condtion  into  a  CFD  code, 
the  inlet  domain  is  extended  axially  past  the  third  stage  rotor.  In  the  axial  direction, 
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the  cross-sectional  area  across  the  first  stage  rotor  is  reduced  linearly  by  relative 
amount  of  area  reduction  calculated  using  Equation  (52).  The  acoustic  disturbance  is 
reflected  off  this  reduced  cross-sectional  area  of  the  domain.  In  addition,  a  constant 
pressure  boundary  condition  is  applied  at  the  exit  of  the  domain.  However,  the 
reflection  off  this  constant  pressure  boundary  was  far  enough  downstream  that  it  did 
not  interfere  with  the  reflected  disturbance  from  the  area  reduction  section  of  the 
domain. 
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4  Governing  Equations 

Both  one-dimensional  and  axisymmetric  two-dimensional  models  were  used  to  numer¬ 
ically  simulate  the  compressor  face  boundary  condition  experiment  at  the  University 
of  Cincinnati.  The  one-dimensional  model  used  the  Euler  equations  in  quasi  one¬ 
dimensional  form  to  calculate  the  unsteady  flow  in  the  inlet.  The  axisymmetric  model 
employed  the  compressible  Favre  averaged  Navier-Stokes  equations  to  investigate  the 
two-dimensional  acoustic  reflection  of  the  axial  compressor. 


4.1  One-Dimensional  Model  Equations 

4.1.1  Euler  Equations 


The  Euler  equations  were  derived  from  the  universal  laws  of  Conservation  of  Mass, 
Conservation  of  Momentum,  and  Conservation  of  Energy  and  governed  the  motion  of 
an  inviscid  fluid.  The  quasi  one-dimensional  Euler  equations,  written  in  conservative 
form,  with  no  body  forces  or  external  heat  addition  were  as  follows. 


<9  (US)  d(F  S)  = 
dt  dx 

where  U  was  the  vector  of  conservative  variables 


W 


U 


P 

pu 

Et 


F  was  the  inviscid  flux  in  the  axial  direction 

pu 

pu2  +p 

[  (Et+p)u 

W  represented  the  source  terms  for  the  equations 


W  = 


0 

vss 

4% 


(53) 


(54) 


(55) 


(56) 


where  S(x,t )  was  the  inlet  cross-sectional  area.  The  source  term  for  the  energy 
equation,  -pff ,  accounted  for  the  work  done  by  the  fluid  on  the  bump  as  the  bump 
deflated. 

4.1.2  Ideal  Gas  Law 

The  three  Euler  equations  contain  three  thermodynamic  variables  p,  p,  Et  as  well  as 
the  axial  velocity  component  u.  Therefore,  another  equation  was  needed  to  close 
the  system.  Air  was  assumed  to  be  a  calorically  perfect  ideal  gas  for  this  work.  Its 
thermodynamic  variables  were  related  by  the  following  equation  of  state. 

P  =  (7  -  1)(E,  ~  \rf)  (57) 

where  Et  was  the  total  energy. 

4.1.3  Geometry  Equations 

For  the  one-dimensional  case,  the  geometry  of  the  experimental  set-up  was  modeled  in 
Equation  (53)  by  the  axial  coordinate  x  and  the  cross-sectional  area  S.  The  constant 
area  annular  duct  containing  the  flexible  bump  was  modeled  using  the  duct  and  bump 
dimensions  published  by  the  University  of  Cincinnati. [9, 10, 11, 12, 16, 17]  In  addition, 
the  constant  area  section  of  the  duct  was  extended  upstream  0.2266  m  to  account 
for  the  bellmouth  entrance  and  0.0319  m  downstream  to  the  mid-plan  of  the  first 
stage  rotor  since  Freund  and  Sajben  had  determined  that  the  acoustic  reflection  in 
the  experiments  originated  at  this  point. [11] 

In  order  to  account  approximately  for  the  pressure  drop  along  the  duct  caused  by 
turbulent  boundary  layer  on  both  the  hub  and  case,  Schlichting’s  flat  plate  turbulent 
boundary  layer  displacement  thickness  model  was  employed. [23]  The  formula  for  the 
displacement  thickness  model  was 
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For  the  calculation  of  the  cross-sectional  area  at  each  axial  location,  the  displacement 
thickness  was  subtracted  from  the  radius  of  the  case  and  added  to  the  radius  of  the 
hub  or  bump.  Figure  18  shows  the  cross-sectional  area  of  the  numerical  domain  prior 
to  the  deflation  of  the  flexible  bump.  The  center  of  the  flexible  bump  was  located  at 
0.8918  m. 


Figure  18:  One-Dimensional  Model  Cross-Sectional  Area 


After  investigating  the  experimental  data  from  the  collapse  of  the  bump,  Slater 
et.  al.  produced  the  following  relation  for  the  height  of  the  center  of  the  collapsing 
bump.  [24] 


h(t)  =  h0 


cos  (if) 


0.8 


(59) 


where  h0  was  the  initial  height  before  deflation  and  T  was  the  duration  time  required 
for  completely  deflating  the  bump.  Slater  also  stated  that  the  bump  collapsed  with  a 
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uniform  displacement  at  each  axial  location.  This  meant  that  the  ends  of  the  bump 
reached  the  hub  first. 

4.2  Axisymmetric  Model  Equations 

4.2.1  Reynolds  and  Favre  Averaging 

In  order  to  account  for  the  effects  of  turbulence  on  the  mean  flow  properties,  the 
Navier-Stokes  equations  were  modified  by  using  Reynolds  and  Favre  averaging. 

For  Reynolds  averaging,  the  variables  can  be  decomposed  into  mean  and  fluc¬ 
tuating  components.  [25]  The  following  is  an  example  of  Reynolds  averaging  of  the 
density. 

p(xt,  t)  =  p{xi)  +  p'{xi ,  t)  (60) 

where  the  mean  quantity  p  is  defined  as 

^  rt+At 

P(*<)=A1™o^y  p(xi,t)dt  (61) 

To  account  for  the  the  density  variations  in  compressible  flows,  Favre  averaging 
was  used  to  simplify  the  averaging  of  the  Navier-Stokes  equations.  For  Favre  averag¬ 
ing,  a  flow  variable  such  as  the  axial  flow  velocity  was  decomposed  into  the  sum  of 
mean  and  fluctuating  parts  as  follows. [25] 

u(xi,  t )  =  u(xi )  +  u"(xi,  t)  (62) 

where  the  Favre  averaged  quantity  u(xi)  was  defined  as 

1  1  7m 

u(xi)  —  -  lim  —  /  pu(xi,t)dt  =  —  (63) 

p  At->oo  /At  p 

The  following  combination  of  Reynolds  and  Favre  averaged  quantities  were  used 

in  the  Favre  averaged  Navier-Stokes  equations. 

P  =  P  +  p'  p  =  p  +  p'  qi  =  cii  +  q'i 
u  =  u  +  u"  V  =  V  +  v"  Tij  =  Tij  +  rb 
ht  =  ht  +  h'l  et  =  et  +  e'l  T  =  T  +  T" 


(64) 
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4.2.2  Navier-Stokes  Equations 

The  Favre  averaged  Navier-Stokes  equations  were  derived  from  the  universal  laws  of 
Conservation  of  Mass,  Conservation  af  Momentum,  and  Conservation  of  Energy  and 
governed  the  motion  of  a  viscous  turbulent  flow.  The  axisymmetric  Favre  averaged 
Navier-Stokes  equations,  written  in  conservative  form,  with  no  body  forces  or  external 
heat  addition  were  as  follows. 


afaU)  9[y(F-Fv)|  d  [y  (G  —  Gv)] 
dt  dx  dy 


where  U  was  the  vector  of  conservative  variables 

P 

u  =  5 


=  W 


F  and  G  were  the  inviscid  fluxes  in  the  axial  and  radial  directions,  respectively 

r  -Pu  i  r  ps  i 


pu 2  +  p 


pv2  +p 


Fv  and  Gv  were  the  viscous  fluxes  in  the  axial  and  radial  directions,  respectively 


0 

0 

T  XX 

0 

< 

II 

TXy 

Txy 

Tyy 

'U'T xx  VT Xy 

i _ 

UT xy  "t"  VT yy  Qy 

W  represented  the  source  terms  for  the  axisymmetric  equations 


W=  _  _ 
p-ree 

0 
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and  y  was  the  radial  distance  from  the  centerline.  The  stress  and  heat  flux  terms  in 
the  above  equations  were  given  by 


^xx  — 
Txy  = 

Tyy  = 
Tee  = 

Qx 

Qy  = 


,  .  ( ndu  dv  v\ 

^+^{2rx-rv-y) 


(m  +  Mr) 


'  du  dv\ 
dy  +  dx) 


2.  .  /  du 

-(y  +  Vr)  [2- 


du  v\ 
dx  y) 
\.  .  f  du  dv 

;(/<  +  «■)  g^  +  2- 

Mr  ^  5T 
dx 


C-n 


M 


+ 


Pr  Pr> 


M  +  Vt\  dT 


Pr  Pri 


dy 


(70) 


where  y  was  the  laminar  viscosity,  yx  was  the  eddy  viscosity,  Pr  was  the  laminar 
Prandtl  number,  Pty  was  the  turbulent  Prandtl  number,  and  cp  was  the  constant 
pressure  specific  heat. 

Physical  domain  grids  usually  have  nonuniform  grid  spacing  in  one  or  more  di¬ 
rections.  Therefore,  it  was  useful  to  transform  the  nonuniform  physical  grid  to  a 
uniformly  distributed  computational  grid.  This  transform  allowed  the  use  of  body- 
fitted  coordinates  which  may  not  be  orthogonal.  The  Navier-Stokes  equations  can 
be  written  for  a  generalized  coordinate  system  (£,77)  using  the  transformation  of  the 
form 


C  =  €{x,y,t) 

v  =  v(x,y,t)  (71) 

T  =  t 


The  differential  expressions  for  this  transformation  to  generalized  coordinates  were 


given  by 
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dt  =  £xdx  +  iydy  +  ttdt 

dr]  =  r]xdx  +  rjydy  +  r]tdt  (72) 

dr  =  dt 


Similarly,  the  differential  expressions  for  the  inverse  transformation  had  the  form 

dx  =  x^dt  +  Xydjj  +  xTdr 

dy  =  yxid£  +  yvdr]  +  yTdr  (73) 

dt  =  dr 


Using  Equations  (72)  and  (73),  the  transformation  metrics  were  defined  as 


tx  — 

ty  =  —XyJ 

Vx  =  -V(J 
Vy  =  X(;J 

where  the  cell  volume  j  was  defined  as 


(74) 


(75) 


J  =  X&r,  -  Xvys 

Applying  the  chain  rule  of  differentiation  to  Equation  (65)  and  using  the  transfor¬ 
mation  metrics,  the  following  axisymmetric  Favre  averaged  Navier-Stokes  equations 
written  in  strong  conservation  law  form  for  the  generalized  coordinate  system  were 
derived. 


d(yfj)  d(yF )  .  d(y 6) 

r\  ”t~  t-\  j~  I  r\ 


dr 


dt 


dy 


W 


(76) 
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where 


U  = 
F  = 
G  = 
W  = 


U 

J 

U^  +  (F-FV)^  +  (G-GV)^ 
J 

U  Vt  +  (F  —  Fv)  rfo  +  (G  —  Gv)  rjy 

J 

W 

J 


4.2.3  Turbulence  Modeling 


(77) 


The  Spalart-Allmaras  one-equation  turbulence  model  was  used  to  determine  eddy 
viscosity. [26]  The  model  included  eight  closure  coefficients  and  three  closure  functions. 
The  eddy  viscosity  was  defined  as 

/j,t  =  pvfv  i  (78) 


The  equation  for  kinematic  eddy  viscosity  was  written  as 
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dv  _  dv 

4”  —  ^1*8^  Ci uiju 
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with  the  following  closure  coefficients  and  functions 
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1/ 


5  —  y/2Qlyfllj  +  ^ 2(pfv2 


(79) 

(80) 
(81) 

(82) 

(83) 

(84) 


The  tensor  fly  =  \(dui/dxy  —  duj/dxi)  was  the  rotation  tensor  and  d  was  distance 
from  the  closest  wall.  This  equation  was  also  transformed  to  generalized  coordinates 
and  solved  in  conjunction  with  the  axisymmetric  Favre  averaged  Navier-Stokes  equa¬ 


tions. 
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4.2.4  Ideal  Gas  Law  and  Sutherland  Law 

The  four  axisymmetric  Navier-Stokes  equations  and  the  Spalart-Allmaras  equation 
contained  seven  thermodynamic  variables  p,p,et,ht,T,  p,  pr  as  well  as  the  velocity 
components  u,  v.  Therefore,  four  more  equation  were  needed  to  close  the  system.  Air 
was  assumed  to  be  a  calorically  perfect  ideal  gas  for  this  work.  Its  thermodynamic 
variables  were  related  by  the  following  equations. 

p  —  pRT 

P*  =  fa2  +  v2)  (85) 

pht  =  pet+p 


where  et  was  the  total  energy  and  ht  was  the  total  enthalpy. 

In  addition,  Sutherland  Law  was  used  to  determine  the  laminar  viscosity  as 


follows.  [27] 


p  =  1.716  x  1(T5 


(_  \  3/2 

T  \  384 

273  )  f  +  111 


(86) 


4.2.5  Geometry  Description 

For  the  axisymmetric  two-dimensional  case,  the  geometry  of  the  experimental  set-up 
was  modeled  in  Equation  (65)  by  the  axial  coordinate  x  and  the  radial  coordinate  y. 
Again,  the  constant  annular  duct  containing  the  flexible  bump  was  modeled  using  the 
dimensions  provided  by  the  University  of  Cincinnati.  In  addition,  the  constant  area 
section  of  the  duct  was  extended  upstream  0.05  m  to  account  for  the  development 
of  the  turbulent  boundary  on  the  bellmouth  entrance  and  0.0319  m  downstream  to 
the  mid-plan  of  the  first  stage  rotor.  Equation  (59)  was  again  used  to  model  the 
deflating  bump  from  the  University  of  Cincinnati  experiment.  The  bump  collapsed 
with  a  uniform  displacement  at  each  axial  location  such  that  the  ends  of  the  bump 
reached  the  hub  first. 
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5  Numerical  Implementation 

In  the  following  sections,  the  numerical  algorithms  employed  for  the  solutions  of  the 
quasi  one-dimensional  Euler  equations  and  the  axisymmetric  Navier-Stokes  equations 
are  presented.  These  include  the  flux  differencing  schemes,  time  integration  schemes, 
and  boundary  conditions  implementations. 


5.1  One-Dimensional  Numerical  Algorithms 


A  cell-centered  finite- volume  approach  was  used  to  solve  the  one-dimensional  model. 
The  domain  interval  [0,  L]  was  divided  into  N  cells  of  length  Ax,  with  cell  centers  %i 
and  cell  faces  ^±1/2  •  The  finite- volume  discretization  of  Equation  (53)  had  the  form 


(SU)"+‘  -  (SU)J 


=  -R(U)" 


R(u)” 


‘^i+1/2^7+1/2  *^-l/2^T-l/2  _  - —  lJl 

Ax  gn  c*n 


where  R(U)  was  the  residual  vector  and  n  represented  the  time  level. 


5.1.1  Roe’s  Flux  Differencing  Scheme 


The  flux  vector  F  in  the  residual  vector  of  Equation  (87)  was  discretized  using  Roe’s 
upwind  method  with  minmod  limited  MUSCL  variable  extrapolation.  [28,  29]  The  cell 
interface  flux  for  Roe’s  method  was  given  by 

Fi+i/2  =  \  [F(ut)  +  F(Ur)  -  |A|(Ur  -  Ut)] (88) 

where  A  =  |jj  was  the  Jacobian  matrix  of  the  Euler  flux  and  the  subscripts  L  and  R 
were  used  to  designate  the  left  and  right  cell  interface  states,  respectively.  The  Roe 
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matrix  A  was  evaluated  using  the  Roe  averaged  variables  given  by 


Ri+ 1/2 
Pi+ 1/2 
tb+1/2 

Hi+ 1/2 


=2 

ai+l/2 


IPR 
PL 
PLRi+ 1/2 
Ri+1/2UR  +  UL 

Ri+ 1/2  +  1 

-Rt+l/2-Hfl  +  Hl 

Ri+ 1/2  +  1 

(7  —  l)(-^i+l/2  —  2^i+l/2) 


(89) 


where  H  was  the  total  enthalpy  defined  by 

7“  IP  2 

The  Roe  vector  in  Equation  (88)  was  evaluated  as  follows. 


|A|(UB-Ut)  =  ^(ft|-9»rrW) 


3= 1 


where  A,-  were  the  eigenvalues  of  the  Roe  matrix 


(90) 


(91) 


Ai  =  u,  X2  =  u  +  a,  A3  =  u  —  a 


(92) 


r(i)  were  the  corresponding  eigenvectors 
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and  dwj  were  the  wave  amplitudes 


dw\  =  Sp  —  dw2  —  Su  +  dw3  =  8u  — 
2, 


5p 

¥ 


pa 


pa 


(93) 


(94) 


A  problem  with  Roe’s  method  is  that  it  can  admit  an  expansion  shock  as  a 
solution  to  the  Euler  equations.  The  easiest  way  to  handle  this  problem  is  to  add  some 
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numerical  dissipation  in  the  calculation  of  the  eigenvalues.  Harten  and  Hyman[30] 
proposed  to  replace  the  absolute  value  of  the  each  eigenvalue  in  Equation  (91)  with 


where  e  was  defined  as  follows. 

e  =  max  0,  (Ai+1/2  -  A*),  (Ai+i  -  Ai+i/2)  (96) 

In  order  to  achieve  a  higher  order  of  accuracy,  variable  extrapolation  of  the 
primitive  variable  values  to  the  left  and  right  cell  interface  states  was  needed.  The 
primitive  variable  vector  was  given  by 

’  P  ' 

V  =  u  (97) 

.  P  . 

The  higher  order  left  and  right  cell  interface  values  of  the  primitive  variables  were 
calculated  using  a  relaxed  minmod  limited  MUSCL  variable  extrapolation.  [31] 

(VL)i+1/2  =  V,  +  i[(l  -  k)FV«  +  (1  +  k)FVJ 
(Vft)j+i/2  =  Vi+1  —  -[(1  +  /c)5-Vj+i  +  (1  —  k)5+Vj+i] 
rvj  =  Vi-Vt-i 
5+Vi  =  Vi+i-Vi 

5~\  =  minmod  (5“V,  5+V,w)  (98) 

(5+V  =  minmod  (5+V,  6~ V,  w) 

minmod(6,  c,  w)  =  ^  (sign(6)  +  sign(c))  min(|fe|,  |wc|) 

1 

K  ~  3 
w  =  1.7 
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The  operators  5+  and  5“  were  the  forward  difference  and  backward  difference  oper¬ 
ators,  respectively.  The  constant  u  was  a  relaxation  factor  to  reduce  the  diffusive 
nature  of  the  minmod  limiter.  For  the  first  order  calculations,  the  left  and  right  cell 
interface  values  were  set  equal  to  the  i  and  i  +  1  cell  centered  values,  respectively. 

5.1.2  Steady  State  Implicit  Integration  Scheme 

For  the  steady  state  time  integration  of  Equation  (87),  the  cross-sectional  area  S  of 
the  annular  duct  varied  only  with  respect  to  the  axial  location  x.  Therefore,  the 
source  term  in  the  energy  equation  was  zero.  In  order  to  reduce  the  time  required  to 
reach  steady  state  convergence,  a  pseudo-transient  continuation  method  was  used  for 
the  implicit  time  integration  of  the  one-dimensional  model. [32,  33,  34]  The  pseudo¬ 
transient  algorithm  had  the  form 

(w  +  £7I)A’”U  =  "R(Ur  (99) 

where  AmU  was  the  conservative  variable  vector  step  size  given  by 

AmU  =  Um+1  -  Um  (100) 

and  m  represented  the  implicit  iteration  level.  This  time  marching  method  sacrifices 
temporal  accuracy  for  rapid  convergence  to  steady  state  by  varying  the  size  of  the  time 
step  A mt.  Kelley’s  implementation  of  the  GMRES  algorithm  was  used  to  calculate 
the  conservative  variable  vector  step  size  for  each  iteration. [35,  36] 

The  residual  vector  on  the  right  hand  side  of  Equation  (99)  was  discretized  using 
Roe’s  flux  difference  scheme.  The  minmod  slope  limiter  in  the  MUSCL  variable 
extrapolation  is  a  source  of  nonsmoothness  for  this  flux  discretization.  Therefore, 
a  smooth  first  order  Lax-Freidrichs  discretization  of  fluxes  was  used  to  construct 
an  approximate  Jacobian  on  the  left  side  of  the  equation.  Lax-Freidrichs  flux 
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discretization  had  the  form 

Pi+i/2  =  2  P^*+i  +  F»  —  |A|j+i/2(Uj+i  —  Uj)]  (101) 

where  |A|  =  |u|  +  a. 

5.1.3  Time  Accurate  Runge-Kutta  Integration  Scheme 


For  the  unsteady  time  integration  of  Equation  (87),  the  new  cross-sectional  area  Sn+1 
over  the  collapsing  bump  was  calculated  using  the  new  bump  height  given  by  Equation 
(59).  In  order  to  insure  the  time  accuracy  of  the  unsteady  inlet  flow  simulation,  an 
explicit,  multi-stage  Runge-Kutta  integration  algorithm  was  chosen. [37]  This  scheme 
provided  higher  order  time  accuracy  by  calculating  the  residual  vector  R  in  Equation 
(87)  at  several  different  values  of  U  in  the  time  interval  At  and  combining  these 
calculations  to  get  Un+1.  The  four-stage  Runge-Kutta  scheme  employed  for  the  higher 
order  calculations  had  the  following  form. 
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uf 

uf 
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At  R(U)f 
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2  sf) 

(SU)» 
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.  R(U)f 

A*  sr1 

(102) 


For  the  first  order  time  accurate  calculations,  an  Euler  explicit  time  integration 
scheme  was  employed. 


5.1.4  Boundary  Conditions 


For  subsonic  inflow,  two  physical  conditions  and  one  numerical  condition  must  be 
specified  at  the  boundary  to  the  domain.  For  each  experimental  run  at  the  Univer- 
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sity  of  Cincinnati,  total  pressure  and  total  temperature  were  recorded.  [11]  Therefore, 
a  constant  total  pressure  and  temperature  were  enforced  physically  at  the  inflow 
boundary.  The  axial  velocity  was  linearly  extrapolated  from  the  domain  at  each  time 
step.  Prom  these  conditions,  the  density  and  pressure  at  the  inflow  boundary  were 
calculated. 

For  a  subsonic  outflow  boundary,  one  physical  condition  and  two  numerical  con¬ 
ditions  must  be  specified.  For  steady  state  time  integration,  a  constant  Mach  num¬ 
ber  boundary  condition  was  enforced.  For  the  unsteady  time  integration,  common 
subsonic  outflow  boundary  conditions  of  constant  Mach  number,  constant  pressure, 
and  constant  axial  velocity  were  employed.  In  addition,  the  Paynter  small  distur¬ 
bance,  Sajben  small  disturbance,  and  new  small  disturbance  boundary  conditions 
were  applied  using  Equations  (10)  and  (48).  The  reflection  of  the  incident  acoustic 
disturbance  created  by  each  of  these  outflow  boundary  conditions  was  compared  with 
the  experimental  data  from  the  University  of  Cincinnati  to  determine  which  one  best 
modeled  the  reflection  characteristics  of  an  axial  compressor. 

5.2  Axisymmetric  Numerical  Algorithms 

A  cell-centered  finite- volume  approach  was  also  used  to  solve  the  axisymmetric  model. 
A  physical  grid  was  developed  for  the  annular  inlet  using  equal  spacing  of  mesh  points 
along  the  axial  direction  and  clustering  of  points  near  the  hub  and  case  for  turbulent 
flow.  This  physical  grid  had  378  x  81  mesh  points.  Figure  19  shows  the  clustering 
of  the  grid  lines  over  the  bump.  Equation  (74)  was  used  to  transform  the  unevenly 
spaced  physical  grid  into  an  evenly  spaced  body-fitted  computational  mesh.  The 
body  surfaces  of  the  hub  and  case  were  t]  =  constant  surfaces.  The  finite- volume 
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Figure  19:  Axisymmetric  Physical  Grid  Over  the  Bump 


discretization  of  Equation  (76)  had  the  form 


vuy  _  (yv\n 

J  )u  \  J  )u 


h3  _ 


At 


=  -R(U) 


n 


(103) 


R(U)J,-  =  (!/F)?+1/2j.  -  (»F)?.1Aj  +  (yO)lj+1/2 
-  -  w?j 

The  indices  (i,j)  represented  the  cell  center  values,  and  the  indices  (i  ±  1/2,  j)  and 
(i,j  ±  1/2)  represented  the  right/left  and  top/bottom  cell  interfaces,  respectively. 


5.2.1  Flux  Representation 

The  flux  definitions  for  F  and  G  given  in  Equation  (77)  contained  inviscid,  viscous, 
and  grid  movement  components  in  the  £  and  r]  directions,  respectively. 

The  inviscid  portion  of  the  flux  vectors  F  and  G  had  the  following  forms.  For 
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clarity,  the  superscripts  denoting  the  Favre-  and  Reynolds-averaged  quantities  will  be 
omitted  for  the  remainder  of  this  dissertation. 


F1  = 
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J 


pU 

pUu  +  ixv 
fMV+iyP 
pUht 

pV 

pVu  +  rjxp 
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(104) 


U  and  V  were  contravariant  velocity  components  in  the  £  and  rj  directions, respectively. 
These  contravariant  velocity  components  were  defined  as 


U  =  £ tu  +  ZyV  (105) 

V  =  rjxu  +  TjyV 


The  inviscid  portions  of  residual  flux  vectors  in  Equation  (103)  were  discretized 
using  Roe’s  upwind  method  with  minmod  limited  MUSCL  variable  extrapolation. [28, 
29]  Roe’s  algorithm  is  an  upwinding  scheme  in  one  dimension  which  is  applied  to 
both  generalized  coordinate  directions  independently.  In  this  section,  Roe’s  scheme 
is  applied  to  the  inviscid  flux  in  f  direction.  The  cell  interface  flux  for  Roe’s  method 
was  given  by 


■£</ 

*  *+1/2  ,j 


Fj(Ul)  +  F7(Ur)  -  I A|(Uh  -  UL) 


j  i+l/2,j 


(106) 


In  order  to  achieve  a  higher  order  of  spatial  accuracy,  variable  extrapolation  of 
the  primitive  variable  vector  to  the  left  and  right  cell  interface  states  was  needed. 
The  primitive  variable  vector  was  given  by 

P 
u 
v 
P 


V  = 


(107) 
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The  higher  order  left  and  right  cell  interface  values  of  the  primitive  variables  were 
calculated  using  Equation  (98)  with  ui  =  1.0. 

The  Roe  vector  in  Equation  (106)  was  given  by 


where 
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The  geometry  and  normal  velocity  terms  in  the  previous  equations  were  defined  as 


|V*| 

,  =  Jv_ 

ivei 

u  =  kxu  +  kyv 


(110) 
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The  Roe  averaged  variables  used  in  Equations  (108)  and  (109)  were  defined  as 


Ri+ 1/2  = 
Pi+ 1/2  = 
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(111) 


Again,  Equation  (95)  was  used  to  handle  this  method’s  problem  of  creating  non¬ 
physical  expansion  shocks. 

For  the  axisymmetric  grid,  a  thin-layer  approximation  was  made  to  the  viscous 
flux  vectors  defined  in  Equations  (70)  and  (77)  by  neglecting  the  derivatives  in  the 
stress  tensor  corresponding  to  the  £  direction.  [31]  In  addition,  the  viscous  flux  in  the  £ 
direction  was  ignored.  Therefore,  the  viscous  portion  of  the  flux  vector  G  in  Equation 
(103)  was  reduced  to  the  following  form. 
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where 
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Central  differences  about  the  cell  interface  were  used  for  the  partial  differentials  in 


Equation  (113).  In  addition,  the  thin  layer  approximation  was  applied  the  stress 
tensor  source  term  in  Equations  (70)  and  (77).  The  thin  layer  source  stress  term  had 


the  form 

1  TL  2. 

-jree  =  jjU*  +  Pr) 

The  grid  movement  portion  of  the  flux  vectors  F  and  G  in  Equation  (77)  con¬ 
tained  the  grid  movement  metrics  ^  and  ‘uj,  respectively.  These  metrics  were  defined 
as  follows 
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5.2.2  Runge-Kutta  Time  Integration 

The  steady  and  unsteady  time  integration  of  Equation  (103)  used  a  two-stage  Runge- 
Kutta  scheme  for  the  higher  order  calculations.  [37]  This  time  integration  scheme  had 
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the  following  form. 
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For  the  first  order  time  accurate  calculations,  an  Euler  explicit  time  integration 
scheme  was  employed. 

5.2.3  Boundary  Conditions 


For  the  axisymmetric  model  inflow,  a  constant  total  temperature  and  mass  flow  rate 
were  enforced  physically  at  the  boundary.  In  addition,  the  radial  velocity  was  set 
equal  to  zero  along  the  boundary,  and  the  pressure  was  linearly  extrapolated  from 
the  interior  of  the  inlet.  From  these  conditions,  the  axial  velocity  and  density  at  the 
inflow  boundary  were  calculated. 

For  steady  state  time  integration,  a  constant  pressure  boundary  condition  was 
enforced  along  the  outflow  boundary.  For  the  unsteady  time  integration,  the  Paynter 
small  disturbance  and  new  small  disturbance  boundary  conditions  were  applied  at  the 
outflow  boundary  using  Equations  (10),  (48),  and  (49).  Again,  the  reflection  of  the  in¬ 
cident  acoustic  disturbance  created  by  each  of  these  outflow  boundary  conditions  was 
compared  with  the  experimental  data  from  the  University  of  Cincinnati  to  determine 
which  one  best  modeled  the  reflection  characteristics  of  an  axial  compressor. 

A  non-slip  wall  boundary  condition  was  applied  at  the  solid  walls  of  the  case  and 
the  hub.  However,  the  normal  non-slip  wall  boundary  condition  was  modified  for  the 
flexible  bump  boundary  to  account  for  the  bump  movement.  This  modified  non-slip 
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boundary  condition  had  the  following  form  along  the  bump. 

V  =  g  (117) 

V  was  the  flow  velocity  at  the  wall,  and  g  =  xti  +  ytj  was  grid  speed  vector  of  the 
collapsing  bump.  This  boundary  condition  provides  the  mechanism  for  the  flow  to 
sense  the  movement  of  the  bump.  In  addition,  the  pressure  and  density  were  extrap¬ 
olated  from  the  interior  domain  to  the  ghost  nodes  along  the  solid  wall  boundaries 
using  a  zero-order  extrapolation. 
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6  Results  and  Discussion 

The  following  sections  discuss  the  results  from  this  current  work.  First,  the  discussion 
will  focus  on  the  compressor  face  boundary  condition  results  for  the  one-dimensional 
model.  These  results  will  be  compared  to  the  experimental  data  from  the  University 
of  Cincinnati.  Then,  the  discussion  will  investigate  the  performance  of  the  boundary 
conditions  for  the  axisymmetric  numerical  model  of  the  experimental  setup. 

6.1  One-Dimensional  Model 

6.1.1  Model  Validation 

Freund  and  Sajben  at  the  University  of  Cincinnati  gathered  experimental  data  on 
the  static  pressure  distribution  across  the  annular  duct  for  an  inlet  Mach  number 
of  0.1 6.  [11]  Table  1  lists  the  flow  and  compressor  conditions  for  this  experimental 
run.  The  listed  ambient  pressure  pamb  and  temperature  Tamb  were  used  for  the  total 
pressure  and  temperature  conditions,  respectively,  at  the  inflow  boundary  for  the 
numerical  one-dimensional  model.  The  reference  pressure  pref  was  measured  at  an 
axial  location  approximately  1.53  meters  in  front  of  the  mean  plane  of  the  first  stage 
rotor.  The  reference  dynamic  pressure  qrej  was  defined  to  be  the  total  pressure  minus 
the  reference  pressure. 


parnb  (kPa) 

98.78 

Tamb  (K) 

300.2 

Ncor  (RPM) 

18452 

pref  (kPa) 

96.92 

Qref  (kPa) 

1.86 

Table  1:  Flow  and  Compressor  Conditions  for  Mq  =  0.16  Experimental  Run 


Figure  20  shows  a  comparison  of  the  static  pressure  distribution  calculated  by 
the  one-dimensional  numerical  model  and  the  distribution  measured  at  the  University 
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of  Cincinnati  for  this  Mach  0.16  experimental  test.  For  this  figure,  x  =  0.0  repre¬ 
sents  the  mid-plane  of  the  first  stage  rotor.  The  flat  plate  turbulent  displacement 
thickness  model  given  in  Equation  (58)  was  responsible  reducing  the  cross-sectional 
area  along  the  inlet.  This  reduction  in  area  caused  the  pressure  drop  across  the  inlet. 
In  Figure  20,  the  pressure  distribution  calculated  by  this  model  agrees  well  with  the 
experimentally  measured  static  pressure  distribution  except  at  the  location  of  the 
struts  in  the  inlet.  Struts  were  located  at  approximately  -0.15,  -1.25,  and  -1.7  meters. 
These  struts  were  not  modeled  in  the  numerical  code.  As  the  results  in  the  next 
sections  will  show,  the  exclusion  of  the  struts  did  not  affect  the  ability  of  this  model 
to  replicate  the  incident  and  reflected  disturbances  travelling  through  the  inlet.  In 
addition,  the  flat  plate  displacement  thickness  model  overpredicted  the  displacement 
thickness  and  subsequently  the  decrease  in  pressure  over  the  flexible  bump  curvature 
as  shown  in  the  figure.  To  account  for  this  overprediction  of  displacement  thickness 
and  pressure  decrease,  the  height  of  the  bump  was  increased  approximately  10%  to 
more  closely  match  the  amplitude  of  the  experimental  incident  expansion  wave  cre¬ 
ated  from  the  deflation  of  the  bump.  Some  separation  right  behind  the  bump  which 
was  not  simulated  in  the  numerical  model  is  shown  in  Figure  20.  This  one-dimensional 
representation  of  the  experimental  setup  at  the  University  of  Cincinnati  provided  a 
good  model  to  test  the  new  compressor  face  boundary  conditions  against  the  data  for 
four  different  experimental  runs. 

6.1.2  Experimental  Test  Run  14A 

First,  the  numerical  one-dimensional  model  was  tested  against  the  experimental  data 
from  the  University  of  Cincinnati  run  14A.[11]  Table  2  lists  the  flow  and  compressor 
conditions  for  this  experimental  test.  For  this  run,  a  Courant-Friedrichs-Lewy  (CFL) 
condition  number  of  0.7  was  employed.  Unless  otherwise  indicated,  all  numerical 
results  presented  for  this  test  case  were  calculated  using  the  higher  order  spatial 


Figure  20:  Axial  Static  Pressure  Distribution  in  the  Inlet  for  M0  =  0.1' 


M0 

0.1637 

Mcf 

0.1867 

Pamb  (kPa) 

Tamb  (K) 

295.7 

2.52 

(RPM) 

17912 

19.8° 

P rotor 

52.0° 

hbump  (m) 

0.0121 

Tbump  (sec) 

0.0008 

Table  2:  Flow  and  Compressor  Conditions  for  Experimental  Run  14A 
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discretiztion  from  limited  MUSCL  variable  extrapolation  and  the  four  stage  Runge- 
Kutta  time  integration  scheme. 

Figure  21  shows  the  steady  state  pressure  and  Mach  number  profiles  for  the  inlet 
prior  to  the  collapse  of  the  bump.  From  a  spectral  analysis  of  their  experimental  data, 
Freund  and  Sajben  concluded  that  the  main  reflection  off  the  compressor  occurred 
at  the  mid-plane  of  the  first  stage  rotor.  [11]  Therefore,  the  domain  extends  to  this 
location  (x  =  2.0349  m)  for  all  steady  state  figures  throughout  the  one-dimensional 
results  section. 


Figure  21:  Steady  State  Pressure  and  Mach  Number  Distributions  for  Run  14A 
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Figure  22  shows  the  creation  of  the  two  expansion  waves  from  the  collapsing 
bump,  the  travel  of  the  waves  through  the  inlet,  and  the  reflection  of  both  waves  off 
the  inflow  and  outflow  boundaries.  The  displacement  of  the  flexible  bump  causes  an 
extraction  of  energy  from  the  flow  and  a  creation  of  a  local  depression  in  pressure. 
After  the  bump  completely  collapses,  this  pressure  depression  splits  into  two  acoustic 
expansion  pulses  which  propagate  in  opposite  directions  through  the  inlet.  After 
the  downstream  travelling  disturbance  interacts  with  the  compressor,  it  is  reflected 
back  upstream  as  an  expansion  pulse  with  a  smaller  amplitude.  For  the  numerical 
results  of  Figure  22,  the  new  small  disturbance  compressor  face  boundary  condition 
was  employed.  At  the  entrance  to  the  inlet,  the  upstream  acoustic  wave  is  reflected 
as  a  compression  wave  since  a  constant  total  pressure  condition  was  enforced  at  the 
inflow  boundary. 

The  accuracy  of  the  downstream  travelling  pulse  in  the  numerical  results  is  of 
primary  importance  to  the  study  of  compressor  face  boundary  conditions.  Figure  23 
shows  the  time  history  of  this  incident  expansion  pulse  as  recorded  by  transducer  sta¬ 
tion  1  (x  =  1.3427  m).  Both  the  first  order  spatial,  Euler  explicit  method  and  higher 
order  spatial,  four  stage  Runge-Kutta  method  created  incident  pulses  that  agreed  well 
with  the  experimental  data  for  the  incident  wave  at  this  station.  Freund  and  Sajben 
stated  that  the  observed  oscillations  behind  the  incident  wave  were  due  to  rebounding 
of  the  flexible  bump  from  the  hub  cage  as  it  deflated. [11]  They  were  unable  to  take 
any  measurements  of  this  effect.  Therefore,  the  rebounding  effect  would  be  nearly 
impossible  to  numerically  simulate  and  was  neglected  for  this  study.  The  amplitude 
of  the  experimental  incident  expansion  pulse  was  -3576  Pa.  The  relative  percent  dif¬ 
ferences  of  the  numerical  amplitudes  compared  with  the  experimental  amplitude  were 
2.29%  and  1.35%  for  the  first  order  and  higher  order  numerical  models,  respectively. 

The  traditional  outflow  conditions  of  constant  pressure,  constant  axial  velocity, 
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Figure  23:  Incident  Acoustic  Expansion  Pulse  at  Station  1  for  Run  14A 

and  constant  Mach  number  were  applied  at  the  compressor  face  boundary.  The  nu¬ 
merical  reflection  characteristics  of  these  boundary  conditions  were  compared  with 
the  experimental  reflection  results  at  transducer  station  4  (x  =  1.9523  m).  Figure  24 
shows  differential  pressure  results  at  station  4  for  the  traditional  boundary  conditions 
and  experiment  for  run  14A.  At  station  4,  the  incident  and  reflected  waves  overlap. 
In  order  to  investigate  the  shape  and  amplitude  of  the  reflected  expansion  pulse,  the 
incident  pulse  measured  at  transducer  1  was  shifted  in  time  and  subtracted  from  the 
pressure  data  collected  at  transducer  4.  Figure  25  shows  the  reflected  acoustic  pulses 
produced  by  the  traditional  boundary  conditions  compared  with  the  axial  compressor 
data  for  this  test  case.  The  constant  pressure  boundary  condition  reflects  a  compres¬ 
sion  wave  instead  of  an  expansion  wave.  The  constant  axial  velocity  and  constant 
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Figure  24:  Comparison  of  Traditional  Boundary  Conditions  vs.  Experimental  Data 
at  Station  4  for  Run  14 A 


Figure  25:  Traditional  Boundary  Conditions’  Acoustic  Wave  Reflections  for  Run  14A 
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Mach  number  conditions  predict  an  expansion  reflection  over  twice  the  amplitude 
measured  during  the  experiment.  In  agreement  with  previous  work,  this  study  also 
showed  that  traditional  outflow  boundary  conditions  do  not  properly  represent  an 
axial  compressor’s  reflection  characteristics. 

Figures  26  and  27  show  the  reflection  characteristics  of  the  Paynter  and  Sajben 
small  disturbance  models  at  transducer  station  4  for  run  14A.  The  largest  amplitude 
of  the  experimental  reflection  pulse  was  measured  as  -1335  Pa  which  corresponded  to 
a  36.8%  reflection  of  amplitude  of  the  incident  pulse  olf  the  axial  compressor  for  this 
test  run.  The  Sajben  small  disturbance  model  predicted  a  reflection  disturbance  that 
was  less  than  half  that  measured  by  the  University  of  Cincinnati.  From  Figure  27, 
the  amplitude  of  the  Paynter  small  disturbance  reflection  pulse  is  -1262  Pa  which  is 
a  5.47%  relative  difference  from  the  experimental  peak  amplitude.  In  addition,  the 
amplitude  of  this  Paynter  reflection  was  35.3%  of  the  numerical  incident  wave  am¬ 
plitude.  Therefore,  the  Paynter  small  disturbance  model  provided  better  agreement 
with  the  initial  and  largest  peak  of  the  experimental  reflection  than  the  Sajben  small 
disturbance  model. 

Due  to  the  proprietary  nature  of  the  blade  configuration  for  the  axial  compressor 
used  at  the  University  of  Cincinnati,  a  few  assumptions  had  to  be  made  with  regard 
to  the  input  values  required  by  the  new  small  disturbance  model.  For  Equation  (26), 
the  ratio  of  cross-sectional  area  in  front  of  the  rotor  to  the  area  behind  was  chosen  to 
be  1.0  because  a  typical  first  stage  rotor  has  little  change  in  area  across  the  blades. 
In  addition,  values  for  the  inlet  guide  vanes  and  first  stage  rotor  exit  flow  angles  were 
needed.  The  known  stagger  angle  for  the  variable  inlet  guide  vanes  was  used  for  ai  in 
Equation  (47)  since  these  vanes  were  uncambered  blades.  The  known  stagger  angle 
for  the  first  stage  rotor  blades  was  used  for  /?4  in  Equation  (47).  The  stagger  angle 
for  the  rotor  is  larger  than  the  exit  blade  angle.  However,  the  exit  angle  of  the  flow 


Figure  26:  Comparison  of  Paynter  and  Sajben  Small  Disturbance  Models  vs.  Exper¬ 
imental  Data  at  Station  4  for  Run  14A 


Figure  27:  Paynter  and  Sajben  Models’  Acoustic  Wave  Reflections  for  Run  14A 
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is  also  larger  than  the  exit  blade  angle.  Therefore,  this  approximation  of  fa  should 
be  fairly  accurate. 

The  numerical  reflection  results  from  the  new  small  disturbance  compressor  face 
boundary  condition  model  were  compared  with  the  experimental  data  from  the  Uni¬ 
versity  of  Cincinnati  run  14A  and  the  numerical  results  from  the  Paynter  small  distur¬ 
bance  boundary  condition.  Figures  28  and  29  show  these  comparisons.  In  Figure  29, 
the  amplitude  of  the  reflected  wave  from  the  new  small  disturbance  model  shows 
excellent  quantitative  agreement  with  the  amplitude  of  the  largest  peak  in  the  exper¬ 
imental  reflected  wave.  Although  there  are  some  differences  in  the  waveform,  the  new 
small  disturbance  reflection  also  shows  good  qualitative  agreement  with  the  experi¬ 
mental  reflection.  Both  the  first  and  higher  order  model  calculations  of  the  reflected 
pulse  from  the  new  small  disturbance  compressor  face  boundary  condition  performed 
better  than  the  Paynter  small  disturbance  results  when  compared  to  the  experimental 
data.  The  higher  order  model  calculations  produced  a  reflected  wave  of  amplitude 
-1315  Pa  which  corresponded  to  a  1.50%  relative  difference  from  the  experimental 
amplitude  of  -1335  Pa.  The  first  order  reflected  wave  had  an  amplitude  of  -1270  Pa. 
In  addition,  the  amplitudes  of  the  first  and  higher  order  reflections  using  the  new 
small  disturbance  boundary  condition  were  35.9%  and  36.8%  of  the  numerical  inci¬ 
dent  wave  amplitudes,  respectively.  Therefore,  the  new  small  disturbance  boundary 
condition  models  the  reflection  characteristics  of  the  axial  compressor  more  closely 
than  Paynter’s  model  for  this  experimental  run. 

The  second  reflection  peak  of  the  experimental  wave  in  Figure  29  was  concluded 
to  be  a  reflection  off  the  mid-plane  of  the  second  stage  rotor  by  Freund  and  Sajben.[ll] 
This  second  reflection  plane  was  not  modeled  in  the  one-dimensional  numerical  code. 
However,  a  numerical  run  was  conducted  to  test  this  conclusion.  From  analyzing 
published  University  of  Cincinnati  figures[ll],  the  axial  location  of  the  mid-plane  of 
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Figure  28:  Comparison  of  New  Small  Disturbance  Boundary  Condition  Model  vs. 
Experimental  Data  at  Station  4  for  Run  14A 


t(sec) 


Figure  29:  New  Small  Disturbance  Model’s  Acoustic  Wave  Reflections  for  Run  14A 
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the  second  stage  rotor  was  estimated  to  be  x  =  2.072  m.  The  numerical  domain  was 
extended  to  the  second  stage  rotor,  and  a  numerical  run  was  conducted  to  investigate 
the  reflection  from  this  plane  using  the  new  small  disturbance  boundary  condition. 
For  this  run,  the  same  reflection  response  coefficient  conditions  used  for  the  first 
stage  rotor  were  employed  for  the  second  stage  rotor  since  no  data  was  available  on 
the  second  stage.  In  addition,  the  presence  of  the  first  stage  rotor  was  ignored  so 
that  the  time  history  of  the  reflection  from  just  the  second  stage  rotor  could  be  ana¬ 
lyzed.  Figure  30  shows  the  time  history  of  the  second  stage  reflection  wave  observed 
at  transducer  station  4.  This  figure  also  shows  that  the  reflection  from  the  second 
stage  rotor  mid-plane  arrives  at  station  4  before  the  second  peak  in  the  experimental 
reflection.  Therefore,  this  second  peak  in  the  experimental  wave  was  reflected  down¬ 
stream  of  the  second  stage  rotor.  From  the  figure,  this  experimental  reflection  took 
place  at  approximately  x  =  2.1314  m.  More  information  on  the  configuration  of  the 
axial  compressor  used  in  the  experiment  would  be  needed  to  better  determine  the 
location  and  cause  of  the  second  peak  in  the  experimental  reflected  wave. 

The  effect  of  changing  the  new  small  disturbance  boundary  condition’s  area  ra¬ 
tio  across  the  rotor  |j-  from  1.0  to  1.1  is  shown  in  Figure  31.  The  fairly  large  10% 
reduction  in  cross-sectional  area  produced  a  19.3%  increase  in  the  amplitude  of  the 
reflection  pulse.  As  stated  previously,  the  cross-sectional  area  typically  changes  very 
little  across  the  first  stage  rotor.  Based  on  the  limited  information  available,  the  as¬ 
sumption  that  the  area  remains  constant  across  the  rotor  will  be  used  for  all  remaining 
numerical  runs. 

The  effect  of  changing  the  exit  flow  angle  /34  for  the  first  stage  rotor  is  shown 
in  Figure  32.  This  5°  change  in  the  flow  angle  produced  an  18.1%  decrease  in  the 
amplitude  of  the  reflected  pulse.  However,  the  stagger  angle  of  the  rotor  blades  is  a 
close  approximation  to  the  flow  exit  angle  since  both  of  these  angles  are  greater  than 
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Figure  30:  Acoustic  Wave  Reflection  from  Second  Stage  Rotor  for  Run  14A 


Figure  31:  Acoustic  Wave  Reflections  for  Multiple  Area  Ratios 
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the  blade  exit  angle.  The  assumption  that  the  exit  flow  angle  is  approximately  equal 
to  the  stagger  angle  will  be  used  for  all  remaining  numerical  runs. 


Figure  32:  Acoustic  Wave  Reflections  for  Multiple  Rotor  Flow  Angles 

Figures  33  and  34  show  the  reflection  characteristics  of  the  simple  area  reduction 
boundary  condition.  The  area  reduction  boundary  condition  for  the  higher  order  nu¬ 
merical  model  produced  a  reflected  expansion  wave  with  an  amplitude  of  -1410  Pa. 
The  first  order  reflection  pulse  had  an  amplitude  of  -1362  Pa.  Although  the  area 
reduction  model  overestimated  the  reflection  characteristics  of  the  axial  compressor 
for  run  14 A,  this  simple  model  still  calculated  the  reflection  amplitude  within  an  ac¬ 
curacy  of  5.5%.  The  numerical  reflection  results  of  run  14A  showed  that  both  the  new 
small  disturbance  model  and  area  reduction  model  compressor  face  boundary  condi¬ 
tions  provided  improved  wave  amplitude  accuracy  over  existing  boundary  conditions 
when  compared  with  the  experimental  data. 
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Figure  33:  Comparison  of  Area  Reduction  Boundary  Condition  Model  vs.  Experi¬ 
mental  Data  at  Station  4  for  Run  14A 


Figure  34:  Area  Reduction  Model’s  Acoustic  Wave  Reflections  for  Run  14A 
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6.1.3  Experimental  Test  Run  010A 

Second,  the  numerical  one-dimensional  model  was  tested  against  the  experimental 
data  from  the  University  of  Cincinnati  run  010A.[11]  Table  3  lists  the  flow  and  com¬ 
pressor  conditions  for  this  experimental  test.  A  CFL  of  0.7  was  also  used  for  these 
numerical  runs.  Again,  all  numerical  results  were  higher  order  calculations  unless 
otherwise  stated.  Figure  35  shows  the  inlet  steady  state  pressure  and  Mach  number 
profiles  for  run  010 A. 


Mo 

0.1637 

Mcf 

0.1869 

Pamb  (kPa) 

100.22 

Tamb  (K) 

294.8 

rnCOT  (kg/m3) 

2.52 

Ncor  (RPM) 

18088 

Fvigv 

24.8° 

r  rotor 

52.0° 

hbump  (m) 

0.0116 

Thump  (SGC) 

0.0008 

Table  3:  Flow  and  Compressor  Conditions  for  Experimental  Run  010A 


Figure  36  shows  the  time  history  of  the  incident  expansion  pulse  as  recorded  by 
transducer  station  1  for  run  010A.  The  amplitude  of  the  experimental  expansion  pulse 
was  -3485  Pa.  Solution  of  the  limited  minmod  variable  extrapolation  Roe  flux  model 
with  a  four  stage  Runge-Kutta  scheme  provided  an  incident  expansion  wave  which 
had  an  amplitude  of  -3484  Pa.  Therefore,  the  higher  order  model  almost  matched  the 
wave  amplitude  exactly.  The  differences  in  the  wave  shapes  between  experimental  and 
numerical  were  due  to  the  bump  collapse  mechanism  model.  The  one-dimensional  nu¬ 
merical  model  accurately  predicted  the  acoustic  expansion  disturbance  which  travels 
downstream  and  interacts  with  the  compressor  face  boundary  condition. 

Figure  37  shows  the  reflected  acoustic  pulses  calculated  by  the  Paynter  and  new 
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Figure  35:  Steady  State  Pressure  and  Mach  Number  Distributions  for  Run  010A 

small  disturbance  boundary  conditions  and  measured  by  the  pressure  transducer  at 
station  4  for  run  010A.  The  largest  amplitude  of  the  experimental  reflection  pulse  was 
measured  as  -1375  Pa  which  corresponded  to  a  39.5%  reflection  of  the  incident  pulse 
amplitude  by  the  axial  compressor.  The  amplitude  of  the  Paynter  model  reflection 
pulse  was  -1229  Pa  which  was  a  10.6%  relative  difference  from  the  experimental  am¬ 
plitude.  In  addition,  the  Paynter  small  disturbance  model  reflected  only  35.3%  of  the 
amplitude  of  the  numerical  incident  wave  from  the  compressor  face  boundary.  Both 
the  first  and  higher  order  new  boundary  condition  model  calculations  of  the  reflected 
pulse  provided  better  agreement  with  the  experimental  data  than  the  Paynter  model 
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Figure  36:  Incident  Acoustic  Expansion  Pulse  at  Station  1  for  Run  010A 

results  as  shown  in  Figure  37.  The  higher  order  new  small  disturbance  boundary 
condition  reflected  an  expansion  wave  of  amplitude  -1355  Pa  which  corresponded  to 
a  1.45%  relative  difference  from  the  experimental  amplitude.  The  first  order  reflected 
wave  had  an  amplitude  of -1311  Pa.  The  first  and  higher  order  models  using  the  new 
small  disturbance  boundary  condition  reflected  37.9%  and  38.9%  of  the  amplitude 
of  the  numerical  incident  waves,  respectively.  Therefore,  the  new  small  disturbance 
boundary  condition  model  again  predicted  the  reflection  characteristics  of  the  axial 
compressor  more  closely  than  the  Paynter  small  disturbance  model. 

Figure  38  shows  the  acoustic  expansion  waves  created  by  the  reflection  of  the 
incident  expansion  pulse  off  the  simple  area  reduction  boundary  condition.  The  area 
reduction  boundary  condition  for  the  higher  order  numerical  model  produced  a  re¬ 
flected  expansion  wave  with  an  amplitude  of  -1288  Pa  which  corresponded  to  a  6.33% 
relative  difference  from  the  experimental  amplitude.  The  first  order  reflection  pulse 
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Figure  37:  Paynter  and  New  Small  Disturbance  Models’  Acoustic  Wave  Reflections 
for  Run  010A 


Figure  38:  Area  Reduction  Model’s  Acoustic  Wave  Reflections  for  Run  010A 
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had  an  amplitude  of  -1243  Pa.  Although  the  area  reduction  model  did  not  perform 
better  than  the  new  small  disturbance  boundary  condition,  this  simple  boundary 
condition  did  predict  the  amplitude  of  the  reflection  pulse  more  accurately  than  the 
Paynter  small  disturbance  boundary  condition.  The  numerical  reflection  results  of 
run  010A  also  showed  that  both  the  new  small  disturbance  model  and  area  reduc¬ 
tion  model  compressor  face  boundary  conditions  provided  improved  wave  amplitude 
accuracy  over  existing  boundary  conditions  when  compared  with  the  experimental 
data. 

6.1.4  Experimental  Test  Run  41A 

Third,  the  numerical  one-dimensional  model  was  tested  against  the  experimental  data 
from  the  University  of  Cincinnati  run  41A.[11]  Table  4  lists  the  flow  and  compressor 
conditions  for  this  experimental  test.  Again,  a  CFL  of  0.7  was  used  for  these  numerical 
runs,  and  all  numerical  results  given  below  were  higher  order  calculations  unless 
otherwise  stated.  Figure  39  shows  the  inlet  steady  state  pressure  and  Mach  number 
profiles  for  run  41A. 


Mo 

0.1724 

Mcf 

0.1968 

Pamb  (kPa) 

98.93 

Tarnb  (K) 

295.6 

mcor  (kg/m3) 

2.65 

Ncot  (RPM) 

18749 

Tv/gv 

24.8° 

P  rotor 

52.0° 

hbump  (m) 

0.0126 

Thump  (S6C) 

0.0008 

Table  4:  Flow  and  Compressor  Conditions  for  Experimental  Run  41 A 


Figure  40  shows  the  time  history  of  the  incident  expansion  pulse  as  recorded  by 
transducer  station  1  for  run  41A.  The  amplitude  of  the  experimental  expansion  pulse 
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Figure  39:  Steady  State  Pressure  and  Mach  Number  Distributions  for  Run  41A 

was  -3719  Pa.  The  collapsing  bump  simulated  by  the  higher  order  numerical  model 
created  an  expansion  disturbance  to  the  flow  with  an  amplitude  of  -3654  Pa.  This 
numerical  expansion  wave  had  a  1.75%  difference  in  amplitude  from  the  experimental 
pulse.  Again,  the  one-dimensional  numerical  model  accurately  predicted  the  incident 
acoustic  expansion  disturbance.  This  is  very  important  since  this  wave  is  reflected  by 
the  compressor  face  boundary  conditions. 

Figure  41  shows  the  time  history  of  the  reflected  acoustic  pulses  calculated  by  the 
Paynter  and  new  small  disturbance  boundary  conditions  at  transducer  station  4  for 
run  41A.  The  largest  amplitude  of  the  experimental  reflection  pulse  was  measured  as 
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Figure  40:  Incident  Acoustic  Expansion  Pulse  at  Station  1  for  Run  41A 

-1342  Pa  which  corresponded  to  a  36.1%  reflection  of  the  incident  pulse  off  the  axial 
compressor  for  this  test  run.  The  Paynter  small  disturbance  boundary  condition  pro¬ 
duced  a  reflected  expansion  pulse  with  an  amplitude  of  -1315  Pa  which  corresponded 
to  a  2.01%  relative  difference  from  the  experimental  amplitude.  In  addition,  the  am¬ 
plitude  of  this  numerical  reflected  pulse  was  35.9%  of  the  amplitude  of  the  numerical 
incident  wave.  The  higher  order  new  small  disturbance  boundary  condition  reflected 
an  expansion  wave  of  amplitude  -1348  Pa  which  corresponded  to  a  0.45%  relative 
difference  from  the  experimental  amplitude.  The  amplitude  of  this  reflection  by  the 
new  small  disturbance  model  was  36.9%  of  the  amplitude  of  the  numerical  incident 
wave.  Again,  the  new  small  disturbance  boundary  condition  provided  an  accurate 
reflection  for  the  compressor  face  boundary. 

Figure  42  shows  the  time  history  of  the  reflected  expansion  pulses  created  by  the 
simple  area  reduction  boundary  condition  at  station  4  for  run  41A.  The  area  reduc- 


80 


Figure  41:  Paynter  and  New  Small  Disturbance  Models’  Acoustic  Wave  Reflections 
for  Run  41 A 


Figure  42:  Area  Reduction  Model’s  Acoustic  Wave  Reflections  for  Run  41A 
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tion  boundary  condition  for  the  higher  order  numerical  model  produced  a  reflected 
expansion  wave  with  an  amplitude  of  -1445  Pa  while  the  first  order  reflection  pulse 
had  an  amplitude  of  -1401  Pa.  Although  the  area  reduction  model  overpredicted 
the  reflection  amplitude,  this  boundary  condition  still  calculated  the  reflection  pulse 
amplitude  within  8%  of  the  experimental  value.  The  numerical  reflection  results  of 
run  41A  showed  the  same  trends  as  the  previous  two  experimental  runs. 

6.1.5  Experimental  Test  Run  40A 

Finally,  the  numerical  one-dimensional  model  was  tested  against  the  experimental 
data  from  the  University  of  Cincinnati  run  40A.[11]  Table  5  lists  the  flow  and  com¬ 
pressor  conditions  for  this  experimental  test.  Like  the  previous  runs,  a  CFL  number 
of  0.7  was  used.  All  numerical  results  for  this  run  were  higher  order  calculations 
unless  otherwise  stated.  Figure  43  shows  the  steady  state  pressure  and  Mach  number 
distributions  for  the  experimental  run  40 A. 


M0 

0.1751 

Mcf 

0.1999 

Pamb  (kPa) 

98.19 

Tarnb  (K) 

292.6 

mCOT  (kg/m3) 

2.69 

Ncor  (RPM) 

18316 

r  vigv 

19.8° 

r  rotor 

52.0° 

hbump  (Hi) 

0.0121 

Thump  (S6C) 

0.0008 

Table  5:  Flow  and  Compressor  Conditions  for  Experimental  Run  40 A 

Figure  44  shows  the  time  history  of  the  incident  expansion  pulse  as  recorded  by 
transducer  station  1  for  run  40A.  The  experimental  expansion  pulse  had  an  amplitude 
of  -3515  Pa.  The  incident  disturbance  created  by  the  higher  order  numerical  model 
had  an  amplitude  of  -3481  Pa  which  corresponded  to  a  0.97%  difference  in  amplitude 
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Figure  43:  Steady  State  Pressure  and  Mach  Number  Distributions  for  Run  40A 

from  the  experimental  pulse.  The  first  order  incident  expansion  wave  had  an  ampli¬ 
tude  of  -3448  Pa.  Like  the  previous  three  cases,  the  higher  order  model  for  run  40A 
predicted  the  amplitude  of  the  incident  pulse  more  accurately  than  the  first  order 
model  as  expected.  The  one-dimensional  modeling  of  the  bump  collapse  produced 
incident  waves  that  agreed  well  with  the  experimentally  measured  waves  for  all  four 
cases. 

Figures  45  and  46  show  the  time  history  of  the  reflected  acoustic  pulses  calculated 
by  the  Paynter  small  disturbance,  new  small  disturbance,  and  area  reduction  bound¬ 
ary  conditions  at  transducer  station  4  for  run  40A.  All  of  these  numerical  boundary 
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Figure  44:  Incident  Acoustic  Expansion  Pulse  at  Station  1  for  Run  40A 

condition  models  overpredicted  the  amplitude  of  the  reflected  expansion  pulse.  The 
amplitude  of  the  experimental  reflection  pulse  was  measured  as  -1114  Pa  which  cor¬ 
responded  to  a  31.7%  reflection  of  the  incident  pulse  by  the  axial  compressor  for  run 
40A.  This  31.7%  reflection  for  the  experimental  case  was  much  lower  than  the  typical 
experimental  reflection  values  of  approximately  36%  observed  by  Freund  and  Sajben. 
A  problem  with  the  data  reduction  for  this  experimental  run  could  have  occurred. 
The  reflected  expansion  pulses  produced  by  the  Paynter  and  new  small  disturbance 
boundary  conditions  had  amplitudes  that  differed  from  the  experimental  amplitude 
by  13%  and  18%,  respectively.  In  addition,  the  area  reduction  boundary  condition 
also  produced  an  expansion  pulse  that  differed  from  the  experimental  pulse  by  18%. 
With  the  exception  of  this  run,  the  new  small  disturbance  and  area  reduction  bound¬ 
ary  conditions  provided  wave  reflections  that  accurately  modeled  the  amplitude  of 
the  reflected  wave  from  the  axial  compressor  within  3%  and  8%,  respectively. 
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Figure  45:  Paynter  and  New  Small  Disturbance  Models’  Acoustic  Wave  Reflections 
for  Run  40A 


Figure  46:  Area  Reduction  Model’s  Acoustic  Wave  Reflections  for  Run  40 A 
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6.2  Axisymmetric  Model 

6.2.1  Model  Validation 

Freund  and  Sajben  measured  a  steady  state  pressure  distribution  along  the  annular 
duct  and  velocity  profiles  both  upstream  and  downstream  of  the  flexible  bump  for  the 
experimental  conditions  given  in  Table  1.  [11]  Figure  47  shows  the  comparison  of  the 
axisymmetric  numerical  pressure  distributions  across  the  inlet  to  the  experimentally 
measured  distribution.  For  all  three  sets  of  data,  the  pressure  was  measured  along 
the  case  of  the  annular  duct.  The  first  order  model  used  an  Euler  explicit  time 
integration  with  a  first  order  spatial  discretization,  and  the  higher  order  model  used 
a  two  stage  Runge-Kutta  time  integration  scheme  with  a  minmod  limited  MUSCL 
variable  extrapolation  for  higher  order  spatial  discretization.  Again,  the  struts  located 
at  approximately -0.15,  -1.25,  and  -1.7  meters  were  not  modeled  in  the  numerical  code. 
The  higher  order  model  showed  good  agreement  in  the  calculation  of  the  pressure 
distribution  along  the  bump  when  compared  with  the  experimental  data.  However, 
both  numerical  models  showed  a  larger  pressure  drop  across  the  bump  than  measured 
experimentally.  Freund  and  Sajben  noticed  that  the  measured  pressure  did  not  fully 
recover  downstream  of  the  bump,  but  they  did  not  invesitgate  the  causes  of  this 
pressure  drop.  [11]  From  analyzing  the  experimental  data,  the  pressure  drop  across 
the  bump  should  be  approximately  70  Pa  if  the  pressure  fully  recovers.  This  pressure 
drop  is  due  to  the  boundary  layer  growth  along  the  bump  portion  of  the  inlet.  The 
University  of  Cincinnati  experimental  data  showed  an  additional  approximately  45 
Pa  pressure  drop  across  the  bump  for  a  total  pressure  drop  of  115  Pa.  The  first 
order  model  calculated  a  very  large  pressure  drop  of  approximately  310  Pa  across  the 
bump,  and  the  higher  order  model  calculated  this  pressure  drop  to  be  approximately 
160  Pa.  Therefore,  the  higher  order  pressure  distribution  had  a  pressure  drop  that 
was  approximately  45  Pa  larger  than  the  experimental  drop.  The  cause  of  this  larger 
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Figure  47:  Axial  Static  Pressure  Distribution  Comparison  Along  Inlet  Case 

pressure  drop  maybe  due  to  the  high  pressure  regions  calculated  by  the  axisymmetric 
numerical  model  at  the  leading  and  trailing  edges  of  the  bump  as  shown  in  Figure  48. 
These  non-physical  regions  are  due  to  excess  diffusion  (dissipation  and  dispersion)  [31] 
produced  by  the  numerical  method  at  low  flow  speeds.  To  reduce  this  excess  numerical 
diffusion,  a  preconditioned  numerical  model  would  need  to  be  employed.  However,  a 
preconditioned  model  would  not  allow  for  time  accurate  tracking  of  the  acoustic  pulses 
generated  by  the  collapse  of  the  flexible  bump.  Since  the  purpose  of  this  research  was 
to  accurately  model  the  reflection  characteristics  of  the  axial  compressor  in  time, 
the  higher  order  non-preconditioned  model  presented  in  Section  5  was  used  for  the 
axisymmetric  results.  In  addition,  the  pressure  distribution  along  the  hub/bump  wall 
had  an  oscillatory  effect  at  the  leading  and  trailing  edges  of  the  bump.  Possible  causes 
for  these  oscillations  are  numerical  dispersion  from  the  higher  order  model  and  the 
thin  layer  approximation  to  the  viscous  flux  terms.  The  effects  of  the  larger  numerical 
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pressure  drop  across  the  bump  and  the  oscillations  in  the  flow  on  the  unsteady  results 


will  be  shown  later  in  this  results  section. 
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Figure  48:  Numerical  Pressure  and  Axial  Velocity  Contour  Plots  over  the  Flexible 
Bump  for  the  Validation  Case 


The  velocity  profiles  for  this  experimental  steady  state  test  case  were  measured 
at  x  =  -0.2218  and  -1.5585  meters.[ll]  Figure  49  shows  the  agreement  of  the  numer¬ 
ical  velocity  profiles  with  the  measured  experimental  profile  upstream  of  the  flexible 
bump.  The  first  and  higher  order  numerical  profiles  overlap  at  this  location.  The 
numerical  model  calculated  the  boundary  layerto  be  slightly  thicker  on  the  hub  and 
too  thin  on  the  case  when  compared  to  the  experimental  data.  However,  the  veloc¬ 
ity  profiles  do  agree  reasonably  well  with  the  experimental  results.  Figure  50  shows 
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the  experimentally  measured  and  numerically  calculated  profiles  downstream  of  the 
bump.  Both  numerical  models  calculated  thicker  boundary  layer  at  the  hub  and  a 
thinner  boundary  layer  at  the  case  than  experiment  for  this  location.  The  higher 
order  velocity  profile  is  approaching  a  fully  developed  profile  for  the  annular  duct. 
The  experimental  velocity  profile  shows  very  little  growth  in  the  hub  boundary  layer 
along  the  inlet.  In  addition,  the  experimental  boundary  layer  on  the  case  occupies 
approximately  60%  of  the  total  annulus  height.  The  boundary  condition  at  the  com¬ 
pressor  face  for  the  steady  state  numerical  calculations  was  constant  static  pressure. 
There  may  be  some  factor  at  the  axial  compressor  that  is  not  being  taken  into  con¬ 
sideration  by  the  model  boundary  conditions  that  contributes  to  this  experimental 
velocity  profile.  However,  no  mention  of  any  additional  factors  was  given  by  Freund 
and  Sajben  other  than  possible  blockage  by  the  struts.  Therefore,  the  axisymmetric 
model  boundary  conditions  were  not  changed. 


Figure  49:  Velocity  Profiles  at  x  =  -1.5585  m  for  the  Validation  Case 
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Figure  50:  Velocity  Profiles  at  x  =  -0.2218  m  for  the  Validation  Case 

The  effectiveness  of  the  Spalart-Allmaras  turbulence  model  was  evaluated.  Fig¬ 
ure  51  shows  a  comparison  of  the  turbulence  model  to  the  Law  of  the  Wall  at  locations 
both  upstream  and  downstream  of  the  flexible  bump.  At  the  upstream  location,  the 
turbulent  boundary  was  still  developing.  The  numerical  velocity  profiles  at  both  lo¬ 
cations  show  good  agreement  with  the  Law  of  the  Wall.  In  addition,  the  cell  center  of 
the  first  cell  off  the  hub  or  case  had  a  y+  value  of  approximately  1.5.  Therefore,  the 
physical  grid  was  partitioned  correctly  for  turbulence  calculations,  and  the  turbulence 
model  performed  properly  for  the  axisymmetric  numerical  code. 

The  compressor  face  boundary  reflection  results  from  axisymmetric  turbulent 
flow  model  were  compared  to  the  experimental  data  gathered  for  four  test  cases. 
The  results  presented  in  the  following  sections  were  calculated  using  the  higher  order 
axisymmetric  model. 
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Figure  51:  Profile  of  Turbulent  Boundary  Layer  at  x  =  -0.2218  and  -1.5585  m 

6.2.2  Experimental  Test  Run  14A 

First,  the  results  from  the  axisymmetric  model  were  compared  to  the  experimental 
data  from  run  14A.  Table  2  lists  the  flow  and  compressor  conditions  for  this  experi¬ 
mental  test.  The  rotor  stagger  angle  listed  in  the  table  is  the  angle  at  the  mean  radius. 
For  all  axisymmetric  model  numerical  runs,  the  stagger  angle  for  the  first  stage  ro¬ 
tor  varies  linearly  from  45°  at  the  hub  to  60°  at  the  case.  Like  the  one-dimensional 
model,  assumptions  had  to  be  made  with  regard  to  the  input  values  required  by  the 
new  small  disturbance  boundary  condition  for  the  axisymmetric  model.  Again,  the 
ratio  of  cross-sectional  area  in  front  of  the  rotor  to  the  area  behind  was  chosen  to 
be  1.0,  the  known  stagger  angle  for  the  variable  inlet  guide  vanes  was  used  for  »i, 
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and  the  known  variation  of  the  stagger  angle  for  the  first  stage  rotor  blades  was  used 
for  /?4.  These  same  assumptions  were  used  for  all  numerical  runs  presented  in  the 
axisymmetric  results  section. 

Figures  52  and  53  show  contour  plots  and  pressure  profiles  of  the  acoustic  wave 
motion  through  the  inlet,  respectively.  The  collapse  of  the  flexible  bump  creates  two 
expansion  pulses,  one  propagating  upstream  and  the  other  downstream  toward  the 
compressor  face.  Then,  both  waves  are  reflected  off  the  inflow  and  outflow  boundaries. 
Both  figures  show  the  oscillatory  nature  of  the  flow  field  between  the  two  expansion 
pulses  at  t  —  1.45  and  2.38  milliseconds.  The  steady  state  pressure  profile  in  Figure  53 
has  a  170  Pa  pressure  drop  across  the  bump.  The  normal  pressure  drop  due  to 
the  growth  of  the  boundary  layer  along  the  bump  should  be  approximately  80  Pa. 
The  additional  90  Pa  pressure  drop  could  be  responsible  for  some  of  the  oscillatory 
behavior  in  the  flow  field  since  the  flow  needs  to  adjust  itself  so  that  its  pressure  profile 
is  continuous  across  the  domain.  However,  a  detailed  look  at  the  flow  field  around  the 
bump  after  it  had  collapsed  showed  that  a  compression  wave  was  formed  along  the  hub 
from  the  downward  movement  of  the  fluid  as  the  bump  collapsed.  This  compression 
wave  merged  with  the  high  pressure  region  at  the  trailing  edge  of  the  bump  and 
propagated  across  the  annulus.  When  the  compression  wave  reached  the  case,  it 
reflected  off  the  case  and  spread  out.  This  compression  wave  then  propagated  back 
across  the  annulus,  reflected  off  the  hub,  and  spread  out.  As  the  compression  wave 
spread  out,  it  produced  oscillations  in  the  flow  that  propagated  both  upstream  and 
downstream  through  the  duct.  The  pressure  profiles  at  y  =  0.099  m  in  Figure  53  show 
that  the  oscillations  were  being  damped  out  over  time.  These  oscillations  do  interact 
with  the  reflected  wave  from  the  compressor  face  boundary  condition.  The  reflected 
wave  at  t  =  3.99  milliseconds  in  both  figures  was  created  by  the  new  small  disturbance 
boundary  condition.  In  Figure  52,  the  incident  expansion  waves  travelling  in  both 
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directions  are  essentially  planar.  This  planar  nature  of  the  flow  validates  the  use  of  a 
one-dimensional  model  for  this  inlet. 

Figure  54  shows  the  time  history  of  the  downstream  travelling  incident  expan¬ 
sion  pulse  as  recorded  by  transducer  station  1.  The  amplitude  of  the  experimental 
incident  expansion  pulse  was  -3576  Pa,  and  the  amplitude  of  the  incident  pulse  for 
the  axisymmetric  model  was  -3295  Pa.  The  relative  percent  difference  of  the  numer¬ 
ical  amplitude  compared  to  the  experimental  amplitude  was  7.85%.  The  oscillatory 
nature  of  the  flow  field  between  the  two  expansion  pulses  probably  contributed  to 
the  reduced  amplitude  of  the  downstream  propagating  numerical  pulse.  Although 
the  numerical  incident  expansion  pulse  had  an  amplitude  less  than  experimental,  this 
pulse  was  still  used  to  test  the  compressor  face  boundary  conditions.  The  percentage 
of  reflected  pulse  amplitude  to  incident  pulse  amplitude  for  the  numerical  results  was 
compared  to  the  percentage  measured  from  the  experimental  results.  As  stated  in 
the  one-dimensional  results,  the  oscillations  in  the  experimental  results  were  thought 
to  be  caused  by  rebounding  of  the  flexible  bump  from  the  hub  cage.  The  oscillations 
observed  in  the  numerical  results  were  similiar  to  the  waves  in  the  experiment,  but 
the  numerical  oscillations  had  a  higher  frequency  and  less  damping  than  the  ones 
measured  in  the  experimental  test  run. 

The  numerical  reflection  results  from  the  new  small  disturbance  compressor  face 
boundary  condition  model  were  compared  with  the  experimental  data  for  the  Univer¬ 
sity  of  Cincinnati  run  14A.  Figure  55  shows  the  comparison  of  the  numerical  results 
with  the  experimental  data  measured  at  transducer  station  4.  The  essential  overlap 
of  the  experimental  results  at  both  the  hub  and  case  illustrate  the  planar  nature  of 
the  acoustic  waves  travelling  through  the  inlet.  The  numerical  results  also  show  the 
same  character.  At  station  4,  the  incident  and  reflected  waves  overlap.  Previously, 
the  incident  pulse  measured  at  transducer  1  was  shifted  in  time  and  subtracted  from 
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Figure  54:  Comparison  of  Incident  Acoustic  Expansion  Pulses  at  Station  1  for  Run 
14A 

the  pressure  data  collected  at  transducer  4  to  investigate  the  amplitude  and  shape  of 
the  reflected  wave.  However,  the  high  frequency  oscillatory  nature  of  the  flow  behind 
the  incident  wave  made  the  calculation  of  the  numerical  reflected  wave  amplitude 
extremely  difficult.  Therefore,  the  pressure  profile  at  t  =  3.99  x  10-3  sec  in  Figure  53 
was  used  to  show  the  reflected  wave  amplitude  and  shape. 

Figure  56  shows  the  reflected  acoustic  pulses  created  by  Paynter  small  distur¬ 
bance,  new  small  disturbance,  and  area  reduction  boundary  conditions  for  run  14A. 
The  damped  oscillations  produced  by  the  oscillatory  nature  of  the  flow  field  between 
the  incident  pulses  interacts  with  the  reflected  pulses.  The  amplitude  of  the  Paynter 
model  reflection  pulse  was  -1112  Pa  which  corresponded  to  a  33.7%  reflection  of  the 
numerical  incident  pulse  amplitude  from  the  compressor  face  boundary  condition. 
The  new  small  disturbance  boundary  condition  produced  a  reflected  wave  of  ampli¬ 
tude  -1161  Pa  which  corresponded  to  a  35.2%  reflection  of  the  numerical  incident  wave 
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Figure  55:  Comparison  of  Pressure  Data  at  Transducer  Station  4  for  Run  14A 


Figure  56:  Reflected  Acoustic  Waves  from  Compressor  Face  Boundary  Conditions  for 
Run  14A 


97 


amplitude.  Finally,  the  area  reduction  reflection  pulse  had  an  amplitude  of  -1264  Pa 
which  was  a  38.4%  reflection.  The  experimental  results  showed  a  36.8%  reflection 
from  the  compressor  face.  All  of  the  boundary  conditions  produced  reflected  pulse 
amplitudes  that  agree  well  qualitatively  with  the  measured  amplitude  of  compressor 
reflection.  Like  the  one-dimensional  results,  the  new  small  disturbance  boundary  con¬ 
dition  performed  better  than  the  Paynter  model  when  compared  to  the  experimental 
amount  of  reflection. 

6.2.3  Experimental  Test  Run  010A 

Second,  the  results  from  the  axisymmetric  model  were  compared  to  the  experimental 
data  from  run  010A.  Table  3  lists  the  flow  and  compressor  conditions  for  this  experi¬ 
mental  test.  Again,  the  rotor  stagger  angle  listed  in  the  table  is  the  angle  at  the  mean 
radius.  Figure  57  shows  steady  state  pressure  profile  along  the  case  for  run  10A.  The 
results  show  a  pressure  drop  of  180  Pa  across  the  bump.  The  normal  pressure  drop 
due  to  the  growth  of  the  boundary  layer  should  be  approximately  80  Pa. 

Figure  58  shows  the  time  history  of  the  incident  expansion  pulse  as  recorded  at 
station  1.  The  amplitude  of  the  experimental  incident  expansion  pulse  was  —3485  Pa. 
The  amplitude  of  the  incident  pulse  for  the  axisymmetric  model  was  —3214  Pa  which 
corresponded  to  a  relative  percent  difference  of  7.78%.  Again,  the  oscillatory  nature 
of  the  flow  probably  contributed  to  this  reduction  in  the  amplitude  of  the  acoustic 
expansion  pulse.  The  amplitude  of  the  numerical  incident  pulse  was  used  to  determine 
the  amount  of  reflection  caused  by  the  compressor  face  boundary  conditions. 

Figure  59  shows  the  reflected  acoustic  pulses  created  by  Paynter  small  distur¬ 
bance,  new  small  disturbance,  and  area  reduction  boundary  conditions  for  run  010A. 
Again,  the  high  frequency  oscillations  in  the  flow  field  interacted  with  the  reflected 
pulses.  The  amplitude  of  the  Paynter  model  reflection  pulse  was  -1080  Pa  which 
corresponded  to  a  33.6%  reflection  of  the  numerical  incident  pulse,  and  the  ampli- 
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tude  of  the  new  small  disturbance  reflected  wave  was  -1186  Pa  which  corresponded 
to  a  36.9%  reflection.  The  area  reduction  boundary  condition  produced  a  reflected 
wave  with  an  amplitude  of  -1226  Pa  which  corresponded  to  a  38.1%  reflection  from 
the  boundary  condition.  The  experimental  results  showed  a  39.5%  reflection  by  the 
compressor.  Again,  the  new  small  disturbance  boundary  condition  produced  a  re¬ 
flected  pulse  that  agreed  better  with  the  experimental  data  than  the  Paynter  small 
disturbance  boundary  condition.  However,  the  area  reduction  boundary  condition 
performed  better  than  both  of  the  small  disturbance  boundary  conditions  for  this 
test  run. 


Figure  59:  Reflected  Acoustic  Waves  from  Compressor  Face  Boundary  Conditions  for 
Run  010A 

6.2.4  Experimental  Test  Run  41A 

Third,  the  results  from  the  axisymmetric  model  were  compared  to  the  experimental 
data  from  run  41A.  Table  4  lists  the  flow  and  compressor  conditions  for  this  experi- 
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mental  test.  Again,  the  rotor  stagger  angle  listed  in  the  table  is  the  angle  at  the  mean 
radius.  Figure  60  shows  steady  state  pressure  profile  along  the  case  for  run  41A.  This 
pressure  profile  shows  a  215  Pa  pressure  drop  across  the  flexible  bump.  The  normal 
pressure  drop  due  to  the  growth  of  the  boundary  layer  should  be  approximately  90 
Pa. 


Figure  60:  Steady  State  Pressure  Profile  for  Run  41A 

Figure  61  shows  the  time  history  of  the  acoustic  expansion  pulse  as  recorded 
at  station  1.  This  incident  pulse  propagates  downstream  and  interacts  with  the 
compressor  face.  The  amplitude  of  the  experimental  pulse  was  —3719  Pa,  and  the 
amplitude  of  the  numerical  pulse  was  —3371  Pa.  The  relative  percent  difference 
between  the  numerical  and  experimental  pulse  amplitudes  was  9.16%.  This  difference 
in  pulse  amplitudes  is  probably  due  to  the  oscillations  in  the  flow  field.  Again,  the 
amplitude  of  the  numerical  pulse  was  used  to  determine  the  percentage  of  reflection. 
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Figure  61:  Comparison  of  Incident  Acoustic  Expansion  Pulses  at  Station  1  for  Run 
41 A 

Figure  62  shows  the  reflected  acoustic  pulses  created  by  Paynter  small  distur¬ 
bance,  new  small  disturbance,  and  area  reduction  boundary  conditions  for  run  41A. 
The  experimental  data  showed  a  36.1%  reflection  from  the  axial  compressor  for  this 
run.  The  new  small  disturbance  compressor  face  boundary  condition  reflected  an 
expansion  pulse  with  an  amplitude  of  -1258  Pa.  This  reflection  was  37.3%  of  the 
incident  numerical  pulse  amplitude.  The  amplitude  of  the  area  reduction  reflection 
pulse  was  -1353  Pa  which  corresponded  to  a  40.1%  reflection,  and  the  amplitude  of 
the  Paynter  model  reflection  pulse  was  -1144  Pa.  This  amplitude  was  a  33.9%  reflec¬ 
tion  of  the  numerical  incident  pulse  amplitude.  The  new  small  disturbance  boundary 
condition  provided  a  reflection  that  agreed  quantitatively  within  2%  of  the  measured 
reflection  characteristics  of  the  compressor. 
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Figure  62:  Reflected  Acoustic  Waves  from  Compressor  Face  Boundary  Conditions  for 
Run  41A 

6.2.5  Experimental  Test  Run  40A 

Finally,  the  results  from  the  axisymmetric  model  were  compared  to  the  experimental 
data  from  run  40 A.  Table  5  lists  the  flow  and  compressor  conditions  for  this  exper¬ 
imental  test.  Again,  the  rotor  stagger  angle  listed  in  the  table  is  the  angle  at  the 
mean  radius.  Figure  63  shows  steady  state  pressure  profile  along  the  case  for  run 
40A.  A  205  Pa  pressure  drop  across  the  bump  was  calculated  for  this  numerical  case. 
Again,  the  normal  pressure  drop  due  to  the  growth  of  the  boundary  layer  should  be 
approximately  90  Pa. 

Figure  64  shows  the  time  history  of  the  acoustic  expansion  pulse  as  calculated  by 
the  transducer  at  station  1.  The  amplitude  of  the  experimental  pulse  was  —3515  Pa. 
The  amplitude  of  the  numerical  incident  pulse  was  —3214  Pa  which  corresponded 
to  a  relative  percent  difference  of  8.56%.  For  all  four  experimental  test  runs,  the 
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Figure  63:  Steady  State  Pressure  Profile  for  Run 
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oscillatory  behavior  of  the  flow  probably  caused  the  amplitude  of  the  numerical  pulse 
to  be  approximately  8%  less  than  the  experimental  pulse  amplitude. 

Figure  65  shows  the  reflected  acoustic  pulses  created  by  the  three  compressor 
face  boundary  conditions  for  run  40A.  The  experimental  data  showed  a  31.7%  re¬ 
flection  from  the  axial  compressor  for  this  run.  Again,  this  reflection  percentage  was 
much  lower  than  that  measured  for  all  other  experimental  runs.  The  new  small  dis¬ 
turbance  compressor  face  boundary  condition  reflected  an  expansion  pulse  with  an 
amplitude  of  -1192  Pa  which  corresponded  to  a  37.1%  reflection  of  the  numerical  inci¬ 
dent  wave.  The  amplitude  of  the  Paynter  model  reflection  pulse  was  -1145  Pa  which 
corresponded  to  a  35.6%.  The  area  reduction  reflection  pulse  had  an  amplitude  of 
-1209  Pa.  Like  the  one-dimensional  results,  all  boundary  conditions  overpredicted  the 
amount  of  reflection  for  this  experimental  test  run.  However,  the  new  small  distur¬ 
bance  boundary  condition  calculated  the  amount  of  reflection  from  the  compressor 
face  within  5.4%  of  the  experimental  amount  for  all  runs  investigated.  If  run  40A 
is  omitted,  the  reflection  from  the  new  small  disturbance  boundary  condition  agreed 
within  2.6%  of  the  experimental  reflection  compared  to  5.9%  for  the  Paynter  small 
disturbance  model.  The  new  small  disturbance  boundary  condition  provided  more 
accurate  reflection  results  for  both  one-dimensional  and  axisymmetric  models  of  the 
University  of  Cincinnati  experiment  than  the  Paynter  boundary  condition  did. 
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Figure  65:  Reflected  Acoustic  Waves  from  Compressor  Face  Boundary  Conditions  for 
Run  40A 
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7  Conclusions 

The  new  small  disturbance  and  area  reduction  compressor  face  boundary  conditions 
were  developed  to  properly  model  the  interaction  of  an  acoustic  pulse  and  an  axial 
compressor.  The  new  small  disturbance  boundary  condition  included  parameters  ne¬ 
glected  in  the  Paynter  small  disturbance  formulation.  This  new  model  accounted  for 
the  inlet  guide  vanes,  cambered  rotor  blades,  and  the  rotational  speed  of  the  com¬ 
pressor.  Therefore,  the  new  small  disturbance  boundary  condition  provided  a  more 
accurate  representation  of  an  axial  compressor.  The  area  reduction  boundary  condi¬ 
tion  used  a  simple  approach.  The  cross-sectional  area  for  the  domain  was  reduced 
across  the  first  stage  rotor  by  an  amount  equal  to  the  amount  of  frontal  area  in  a  blade 
passage  seen  as  closed  to  the  entering  acoustic  wave.  Both  of  these  boundary  condi¬ 
tions  were  applied  in  one-dimensional  and  axisymmetric  models  of  the  experiment  at 
the  University  of  Cincinnati.  Their  acoustic  wave  reflection  results  were  compared  to 
the  experimental  data  measured  for  four  test  cases  and  the  Paynter  small  disturbance 
boundary  condition  results. 

The  one-dimensional  results  showed  that  the  new  small  disturbance  boundary 
condition  accurately  modeled  the  reflection  characteristics  of  the  axial  compressor 
within  6%  for  all  test  cases.  In  addition,  this  new  boundary  condition  provided  im¬ 
proved  reflection  predictions  compared  to  the  Paynter  model.  The  area  reduction 
boundary  condition  did  not  perform  as  well  as  the  new  small  disturbance  bound¬ 
ary  condition.  However,  this  simple  boundary  condition  did  provide  reflections  that 
agreed  well  with  the  experimental  reflection  characteristics  of  the  compressor. 

The  numerical  effects  present  in  the  axisymmetric  model  caused  some  problems 
in  the  calculation  of  the  initial  steady  flow  conditions  and  the  final  evaluation  of  the 
compressor  face  boundary  conditions.  However,  the  reflection  results  showed  that 
the  new  small  disturbance  boundary  condition  still  provided  the  best  agreement  with 
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experimental  data  from  the  University  of  Cincinnati.  Therefore,  the  new  small  dis¬ 
turbance  compressor  face  boundary  condition  was  shown  to  be  the  best  alternative 
for  modeling  the  interaction  between  an  acoustic  disturbance  and  an  axial  compres¬ 
sor.  For  the  axisymmetric  model,  the  area  reduction  boundary  condition  produced 
reflections  that  agreed  well  with  the  experimental  data.  The  easy  of  implementation 
and  good  agreement  with  the  experimental  data  makes  the  area  reduction  model  a 
viable  compressor  face  boundary  condition  alternative  for  design  studies. 

The  one-dimensional  and  axisymmetric  models  provided  very  similar  results  for 
the  reflection  characteristics  of  the  University  of  Cincinnati  axial  compressor.  In 
addition,  the  expansion  pulses  created  by  the  collapsing  bump  for  the  axisymmetric 
model  were  essentially  planar.  Therefore,  the  one-dimensional  model  of  the  flow 
field  through  this  inlet  is  a  valid  and  computationally  less  expensive  alternative  for 
investigating  the  reflection  characteristics  of  the  compressor  face  boundary. 

For  future  research,  both  the  new  small  disturbance  and  area  reduction  boundary 
conditions  need  to  be  implemented  in  inlet  CFD  codes  and  tested  to  see  how  well  they 
can  predict  inlet  unstart.  Based  on  their  performance  compared  to  the  University  of 
Cincinnati  experimental  data,  both  of  these  boundary  conditions  should  provide  more 
accurate  predictions  than  the  currently  imposed  compressor  face  boundary  conditions. 
In  addition,  more  investigation  on  how  to  eliminate  the  numerical  diffusion  in  the 
axisymmetric  flow  model  is  needed.  Possible  solutions  are  development  of  a  time 
accurate  preconditioned  numerical  model  and  use  of  all  viscous  flux  terms  instead  of 
the  thin  layer  approximation  terms. 
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